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ABSTRACT 

We present a study of the natal core harboring the class protostar GF 9-2 
(bolometric luminosity ^ 0.3 Lq, bolometric temperature <^ 20 K; Wiesmeyer et 
al.) in the filamentary dark cloud GF9 (distance = 200 pc) using the Nobeyama 
45 m, CSO 10.4 m telescopes, and the OVRO mm-array. GF9-2 stands unique in 
the sense that it shows H2O maser emission (Furuya et al.), a clear signpost of 
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protostar formation, whereas it does not have a high-velocity large-scale molec- 
ular outflow evidenced by our deep search for ^^CO wing emission. These facts 
indicate that GF9-2 core is early enough after star formation so that it still 
retains some information of initial conditions for collapse. Our 350 fim dust con- 
tinuum emission image revealed the presence of a protostellar envelope with an 
extent of ~ 5400 AU in the center of a molecular core ~ 0.08 pc in size. The mass 
of the envelope is ~ 0.6 Mq from the 350 //m flux density, while LTE mass of the 
core is ~ 3 Mq from N2H+, E.^^CO+, CCS, and NH3 line observations. Combin- 
ing visibility data from the OVRO mm-array and the 45 m telescope, we found 
that the core has a radial density proflle of p{r) oc for 0.003 ^ r/pc ^ 0.08 
region. Molecular line data analysis revealed that the velocity width of the core 
gas increases inward, while the outermost region maintains a velocity dispersion 
of a few times of the ambient sound speed. The broadened velocity width can 
be interpreted as infall. Thus, the collapse in GF9-2 is likely to be described 
by an extension of the Larson-Penston solution for the period after formation of 
a central star. We derived the current mass accretion rate of ~ 3 x 10~^ Mq 
year~^ from infall velocity of ~ 0.3 km at r ~ 7000 AU. Furthermore, we 
found evidence that a protobinary is being formed at the core center. All results 
suggest that GF 9-2 core has been undergoing gravitational collapse for ^ 5000 
years since the formation of central protostar (s), and that the unstable state 
initiated the collapse ~ 2 x 10^ years (the free-fall time) ago. 

Subject headings: ISM: clouds — ISM: evolution — ISM: individual (GF9-2, 
L 1082) — ISM: molecules — stars: formation — stars: prc-main sequence 

1. Introduction 

Stars form inside dense (~ 10^ cm^^) and compact (~ 0.1 pc) molecular cores (e.g., 
Myers & Benson 1983; Mizuno et al. 1994). One of the roles of such cores is to provide 
material for newly forming stars through dynamical accretion caused by gravitational collapse 
of the cores. For isolated low-mass star formation, two extreme paradigms for collapse of 
a core have been proposed. One is a slow, quasi-static model, while the other extreme is 
the dynamic, turbulent one. The former is represented by the Shu's similarity solution (Shu 
1977), and the latter by the Larson- Penston's solution (Larson 1969; Penston 1969) which 
was extended for the period after a protostar formation by Hunter (1977) and Whitworth & 
Summers (1985). Hereafter, we refer to the latter extended solution as the Larson- Penston's 
solution. The two solutions are discriminated by radial density and velocity proflles of the 
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core ait — 0. The Shu's solution (1977) sets the density profile of p(r) = and is static 

(v{r) = 0) at i = where Cg and G are the isothermal sound speed and the gravitational 
constant, respectively. At t = 0, the collapse of the core begins at the center, which produces 
an expansion or rarefaction wave that propagates outward with the local sound speed. This 
model is referred to as the "inside-out" collapse (Shu 1977; Shu, Adams & Lizano 1987). On 

2 

the other hand, the Larson-Penston's solution approaches the profile of p(r) ~ 4.4 
and in infaUing state {v{r) pa —3.28 Cg) when a protostar is about to be born {t — > —0, Hunter 
1977; Shu 1977; Foster & Chevaher 1993). Thus the Larson- Penston solution, referred to as 
the "runaway" collapse, predicts a 48 times higher mass infall rate than the Shu's. At t > 
0, the both solutions have a free-fall region of p{r) oc r~^/^ which expands with sound speed 
for the Shu solution and with supersonic velocity for the Larson-Penston that accompanies 
no rarefaction wave. 

Molecular line and dust continuum observations have identified a number of gravita- 
tionally bound starless cores (e.g., Mizuno et al. 1994), referred to as "preprotostellar" 
cores (Ward- Thompson et al. 1994). Thermal dust continuum emission imaging at mm and 
submm wavelengths revealed that radial density profiles of the preprotostellar cores have 
much flatter profiles than p(r) oc at radii smaller than a few thousand AU (e.g., Ward- 
Thompson et al. 1994; Andre et al. 1996). Molecular line studies measuring velocity field of 
the cores are also powerful tool to verify collapse models (e.g., Zhou 1992). For this purpose, 
it is the most important to search for a star-forming core which retains the initial conditions 
for the gravitational collapse, namely, a star-forming core at an extremely early evolutionary 
phase before dispersal by a large-scale molecular outflow. Furuya et al (1999, 2000) proposed 
that the intermediate-mass class source S106FIR is such an ideal core, because the source 
is characterized by (1) a cold SED (Richer et al. 1993), (2) H2O maser emission, i.e., strong 
evidence for the presence of a protostar, and (3) absence of a large-scale high-velocity CO 
outflow. 

We have successfully identified a low-mass counterpart, GF9-2 by performing multi- 
epoch H2O maser survey (Furuya et al. 2001, 2003) and ^^CO outfiow mapping survey 
toward the selected low-mass protostars with the Nobeyama 45 m, SEST 15 m, and CSO 
10.4m telescopes (R. S. Furuya et al. in preparation). In our sample, GF9-2 is the only 
protostar satisfying the above characteristics of (l)-(3). This protostar is located in the 
east-west filamentary dark cloud GF 9 where several YSOs are located every ~ 0.75 pc in 
the plane of the sky (Schneider & Elmegreen 1979; Wiesemeyer 1997; the filament is cross- 
identified as Barnard 150 and Lyndsl082. The distance to GF9-2 is estimated to be 200 
pc by Wiesemeyer (1997). Wiesemeyer et al. (1997, 1999) identified GF9-2 as class and 
obtained bolometric temperature (Tboi) of < 20 K and bolometric luminosity (Lboi) of O.3L0. 
We believe that GF9-2 offers unique opportunity to investigate the initial conditions for 
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gravitational collapse of a star-forming core. This paper presents a detailed observational 
study of GF9-2. In particular, we combined mm-interferometer and single-dish telescope 
visibility data to obtain high fidelity images of the core. 

2. Observations 

Toward GF9-2, we observed dust continuum emission with the Caltech Submillimeter 
Observatory ^ (CSO) 10.4 m telescope (§2.1) and molecular line emission with the Nobeyama 
Radio Observatory^ (NRO) 45 m telescope (§2.2). We also carried out a deep search for 
molecular outflow using the 10.4 m and 45 m telescopes and used cm radio continuum data 
from the NRAO'^ Very Large Array (VLA) archive data base to search for thermal radio jet 
(§2.4). The parameters for single-dish and interferometer observations are summarized in 
Tables 1 and 2, respectively. 

2.1. Caltech Submillimeter Observatory 10.4 m Telescope 

2.1.1. Mapping of^^CO (3-2) Line Emission 

The ^^CO (3-2) line observations were done using the SIS receiver with position-switching 
mode in 2002 December. We used the 50 MHz bandwidth acousto-optical spectrometer 
(AOS). During the observations, the zenith opacities at 225 GHz (T225 ghz) ranged from 
0.037 to 0.12. The telescope pointing was checked by observing the ^^CO (3-2) emission 
toward CRL 2688 every 1.5 hours, and the accuracy was better than 5". The observing grid 
was centered on the position of the 3 mm continuum source reported in Wiesemeyer et al. 
(1999) (R.A.J2000 = 20''51'"29'.'5, Dec.jsooo = 60°18'38'.'0). In the data reduction, we removed 
spectra obtained with elevation angle of ^ 30° or with Tgys > 800 K in single sideband (SSB). 



^The Caltech Submillimeter Observatory is operated by the California Institute of Technology under 
funding from the US National Science Foundation (AST 02-29008). 

^Nobeyama Radio Observatory is a branch of the National Astronomical Observatory, operated by the 
Ministry of Education, Culture, Sports, Science and Technology, Japan. 

^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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2.1.2. Imaging of SSOjim Continuum Emission with SEARCH 

The SHARC II camera (Dowell et al. 2003) was used to obtain a 350 fim continuum 
emission map in the box and Lissajous scan modes (Kovacs 2006) in 2004 June. We used the 
Dish Surface Optimization System (DSOS) which actively corrects gravitational deformation 
of the main reflector of the telescope (Leong 2005). During our observations, T225 ghz ranged 
from 0.045 to 0.08. The telescope focusing and pointing were checked every an hour by 
observing Neptune and Uranus. The measured beam size in HPBW was (8'/9 ± 0'.'2) x 
(S'.'O ± 0'.'4) at P.A.= 93°. We estimated the positional accuracy of ^ 4" from maps of the 
pointing sources. For flux calibrations, we used the pointing data of Neptune that satisfied 
the condition of T225GHZ X airmass < 0.1. We obtained a mean peak voltage (Vpcak) of 
8.25 ± 1.5 /iVolt toward Neptune over a pixel aperture (Apixci = 4'.'9 x 4'.'8) of the detector. 
The measured voltage was converted into source intensity as oc ^ where Si, denotes 
flux density per beam and beam solid angle. Since the brightness temperature of 63 K 
and the angular diameter of 2.6" for Neptune led to S„ of 98.0 Jy at 350 //m, we obtained a 
conversion factor of 3.2x10"^ Jy beam~^ Volt~^. 

2.2. Nobeyama 45 m Telescope 

2.2.1. Observations of^^CO (1-0), H^^CO+ and HC^^O+ (1-0), and C3H2 2i2 - Iqi Lmes 

We carried out simultaneous observations of ^^CO (1-0) and H^^CO+ (1-0) fines using 
two SIS receivers (SlOO and S80) witfi SSB filters in dual polarization mode in 2000 April. 
For the backend, we used the 8 arrays of AOS. All of the 45 m telescope observations were 
made in position switching mode, using an off-position of (Aa, AS) = {—8'.0, +10^0) to the 
3 mm source. The pointing accuracy was better than 3". The daily variation of the receiver 
gains was estimated to be less than 10% in peak-to-peak, leading the corrected absolute fiux 
uncertainty of 15%. In 2005 April, we made the simultaneous observations of the H^^CO+ 
and HC^^O"*" (1-0) lines to measure optical depth of the H^^CO+ line toward the 3 mm 
continuum source. At the same time, we also observed the ortho-C3H2 2i2 — Iqi line. 

2.2.2. Mapping Observations of^^CO (1-0), H^^CO+ (1-0) and SiO (2-1) Lines with 

BEARS 

The mapping observations of the ^^CO (1-0), Y{^^CO+ (1-0) and SiO (2-1) ^; = 
lines were carried out with the 25 BEam Array Receiver System (BEARS) in 2001 January 
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and April. For backend, we used digital spectrometers. We adopted correction factors 
provided by the observatory for relative gain differences among the 25 beams. The ^^CO 
and H^'^^CO^ line data taken with the S100/S80 receivers with the SSB filters (§2.2.1) were 
used to determine the antenna temperature (T^) scale for the BEARS which is not equipped 
with SSB filters. The telescope pointing accuracy was better than 3", and the uncertainty 
of temperature scale calibration was better than 20%. 

2.2.3. Observations of N2H+ (1-0), CCS As - 82 and CC^^S - 82 Lines 

In 2003 December, we made the simultaneous N2H+ (1-0) and CCS 43 — 82 line obser- 
vations with a grid spacing of 40", full-beam sampling for the CCS line. Subsequently, we 
carried out deeper integration of the N2H'^ line with a spacing of 20", a full-beam sampling 
for the line. In addition, we observed CC^^S 43 — 82 line toward the 8 mm continuum source 
to estimate optical depth of the CCS line. The pointing accuracy and the uncertainty of flux 
cahbration were the same as in §2.2.2. 

2.2.4. Observations of NH^ Inversion Lines at 23 GHz 

We carried out mapping observations of the four inversion lines of NH3 molecules at 
28-24 GHz using the HEMT receiver (H20) with a dual circular polarization mode in 2001 
April. We used the eight AOS arrays to receive both the polarizations of the 4 transitions 
and added the dual polarization data for each transition to improve S/N ratio. The telescope 
pointing accuracy was better than 5", and the uncertainty of temperature calibration was 
better than 20%. 



2.3. Owens Valley Radio Observatory Millimeter Array 

The aperture synthesis observations at 8 mm were carried out using the six-element 
OVRO mm array in the period from 2008 September to 2004 May with the 4 array con- 
figurations; 8 full-tracks at the C configuration, 5 at L, 2 at H, and 2 partial tracks at E 
(Table 2). The resultant projected baseline length ranged from 12.8 m to 224 m; our OVRO 
observations are insensitive to structures extended more than 54", corresponding to 0.052 pc 
at d = 200 pc. We used the 8 mm continuum source position reported in Wiesemeyer et al. 
(1999) as the phase tracking center. We observed the N2H"'' (1-0) line in the upper sideband 
and the H^^CO^ (1^0) liiis in the lower sideband. The velocity resolutions were 0.402 and 
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0.108 km for the N2H+ and R^^CO+ lines, respectively. 

We used 3C273, 3C 454.3 and 3C84 as passband calibrators and J2038+513 as a phase 
and gain calibrator. The flux densities of J2038+513 were determined from observations of 
Uranus and Neptune. The measured flux densities were fairly stable in the range of 1.37 - 
1.43 Jy at 86.7 GHz and of 1.26 - 1.38 Jy at 93.2 GHz. The uncertainty in our flux calibration 

was estimated to be 10%. To minimize atmospheric phase fluctuation errors, we removed all 
the GF 9-2 data taken in the switching cycles when the coherency of the calibrator was less 
than 70% . The data were calibrated and edited using the MMA and MIRIAD packages. 
For the continuum data, we merged the visibilities in both the sidebands to construct an 
image; the final center frequency is 89.964 GHz. Images of both the line and continuum were 
constructed using the AIPS package. 

Since an interferometer does not have sensitivity to structures extended more than the 
largest fringe-spacing given by the shortest projected baseline length, we used the 45 m 
telescope maps to recover filtered-out emission in our OVRO mm-array observations (Table 
3; Appendix A). After combining interferometer and single-dish visibihties, we GLEANed 
maps using the task IMAGR in the AIPS package. The parameters of the final maps are 
given in Table 3. 

2.4. Very Large Array 

We used the 8435 MHz continuum emission image of GF9-2 from the VLA archive 
database. The data were originally taken by Wiesemeyer et al. in 1997 July with the G- 
array configuration (Table 2). The fiux densities of cahbrators of J2022-I-616 and J2055-I-613 
were determined from observations of the standard calibration sources of 3 G286 and 3 G48, 
and the uncertainty in the final fiux calibration was better than 20%. 

3. Results 
3.1. Continuum Emission 

3.1.1. Non-Detection of Free- Free Emission at 8.4 GHz 

No radio continuum emission from ionized gas is detected toward GF9-2 through our 
VLA G-configuration search at 8.4 GHz (A = 3.6 cm). The attained 3a sensitivity is 0.11 
mJy beam~^, leading to the upper limit of the 3.6 cm continuum luminosity of 0.005 mJy 
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kpc^ within the beam. This upper hmit is lower than the luminosity of the least luminous 
thermal cm continuum emission (see Figure 3 of Cabrit & Bertout 1992; Furuya et al. 2003). 
Since cm continuum emission toward protostars is usually attributed to free-free emission 
from ionized gas produced by strong shocks of protostellar jets, the lack of the cm emission 
in GF9-2 suggests that a jet, if exists, is not so powerful to ionize ambient gas. 



3.1.2. Thermal Dust Emission at 3 mm and 350 jim 

In Figure 1 , we present a comparison between the continuum emission maps at 350 //m 
taken with the SHARC II camera and at 3 mm with the OVRO mm-array. The maps show 
that a 3 mm point source is embedded in the extended 350 jim emission. The peak position 
of the 3 mm emission is R.A.= 20^51"^29."827, Dccl.= 60° 18' 38'.'06 in J2000. The 350 /im 
emission has a size of 33'.'0 x 22'.'5, corresponding to 6600 x4500 AU which is comparable 
to the typical size of a protostellar envelope (< 6000 AU; Mundy, Looney & Welch 2000). 
It is most likely that the 350 //m emission traces a protostellar envelope, whereas the 3 mm 
emission is more compact, perhaps a protoplanetary disk. We assume that a protostar is 
located at the peak of the 3 mm emission. By integrating the emission inside the 3(7 contours, 
we obtain the total flux densities (Table 4). Here we consider that the 3 mm emission has 
negligible contribution from optically thin free-fee emission because our non-detection of the 
8.4 GHz continuum emission gives a 3(7 upper limit of 0.090 mJy at 3 mm for the free-free 
emission by extrapolating with ly oc 

Assuming that the dust continuum emission is optically thin, we estimate the envelope 
mass (Mcnv) from the 350 ^m. total flux density. We used an equation of Mcnv = ^ ^"'(Td) 
where /tj^ is the dust mass opacity coefficient, the dust temperature, and i?i^(Td) the 
Planck function. Assuming that dust and gas are weU coupled, we adopt of 9.5 K, a 
mean excitation temperature (Tgx) obtained from the N2H+ (1-0) data (described in §4.1.1). 
The value of at 350 //m is estimated by the usual form of — Ki,o {p/pq)^ . Here we 
used two sets of the representative Kq and (5 values reported in preprotostellar cores (case 
I) and in envelopes around class O/I protostars (case II), because kq is known to change 
with evolutionary class. Recall that the source is very likely to be in an evolutionary phase 
between a preprotostellar core and a class protostar. Given Kq of 0.005 cm^ g^^ at Vq = 231 
GHz (Preibish et al. 1993, Andre, Ward-Thompson, & Motte 1996) and (3 of 2.0 (Beckwith, 
Henning & Nakagawa 2000, and references therein) for case I, we have n of 0.069 cm^ g~^ at 
850 GHz. On the other hand, given kq of 0.01 cm^ g~^ at 231 GHz (Osscnkopf & Henning 
1994) and (3 of 1.0 (Beckwith et al. 2000) for case II, we have k of 0.037 cm^ g-^ at 850 GHz. 
Finally, we obtain the total masses of 0.4 Mq for case I and 0.8 Mq for case II. The inferred 
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mass range is comparable to the least massive envelope masses in class sources known to 
date, i.e., M^nv ^ 0.5 Mq for HH25MMS-MM and L1157MM (Andre, Ward- Thompson & 
Barsony 1993, hereafter AWB93). Using the source size (Qg) of 957 arcsec^ inside the 3a 
level contour, and M^nv, a mean column density of H2 molecules, A^(H2)g//, is estimated to 
be 4.2 X 10^^ and 7.8 x 10^^ cm~^ for cases 1 and 11, respectively. These are converted into a 
mean optical depth of ~ 7 x 10"'^ for both the cases, justifying our assumption of optically 
thin emission. 



3.2. Molecular Line Emission 

3.2.1. Spectral Profile of the Molecular Line Emission toward the Center 

In Figure 2, we present the spectra of ^^CO (1-0), (3-2), N2H+ (1-0) Fi,F=0,l-l,2, 
HCi«0+ (1-0), C3H2 2i2 - loi, CCS, CC^^S A3 - Sa, and SiO (2-1) ^; = lines 
toward the 3 mm continuum source. Figure 3 represents the N2H+ (1-0) spectrum to show 
all the hyperfine components. We summarize the parameters of these spectra in Table 5. 
Intrinsic velocity width (Afint) in the table indicates the velocity width corrected for line 
broadening (where appropriate) due to instrumental velocity resolution and optical depth 
of the line (Phillips et al. 1979). The ^^CO (1-0) spectrum shows a broad line width of 
~ 2 km s~^ in full- width at half maximum (FWHM). The ^^CO (3-2) spectrum displays a 
blueshifted wing emission up to Vlsr — —6 km s~^ and a less prominent red wing up to ~ 
-1-0.5 km s~^. The spectrum has a central dip at Vlsr = —2.8 km s~^, while the high density 
gas tracers of N2H+, H^^CO+, C3H2, and CCS show their peaks around the dip of the ^^CO 
(3-2) spectrum. Such a coincidence indicates that the two peaks in the ^^CO (3-2) spectrum 
arise from self-absorption, and not from two velocity components along the line-of-sight. At 
^^CO (1-0), the self-absorption is broad; the peak in the redshiftcd side is almost missing, 
and only the blueshifted peak is seen. Here, we define systemic velocity (K;ys) of the cloud 
core as —2.5 km s~^, the peak velocity of the CCS 43 — 32 line, because the line has a single 
sharp peak and is optically thin (described in §4.1.4). In addition, the H^^CO^ and N2H"'" 
lines show similar line profiles in the sense that they have stronger blue peaks, a typical 
signature of infall (e.g.. Walker et al. 1986). The C3H2 spectrum shows a fiat-top shape. 
Lastly, we point out that no molecular line emission is seen around Vlsr — 5.7 km s~^ where 
we detected the weak H2O masers (Furuya et al. 2003). This may suggest that the masers 
are excited in shocked region impacted by a possible protostellar jet because such jet-origin 
water masers in low-mass YSOs usually do not appear around their cloud velocities (Terebey, 
Vogel & Myers 1992; Chernin 1995; Claussen et al. 1998; Moscadelh et al. 2006). 
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3.2.2. Single-Dish Maps of High Density Gas 

Figure 4 illustrates total integrated intensity maps of the high density gas tracers of 
(1-0), CCS 4,3 - 32, N2H+ (1-0), and NH3 hues taken with the 45m telescope. 
We integrated all of the hyperfine components for the nitrogen-bearing molecules having 
hyperfine structures. The four lines show similar overall structures. It is very likely that all 
of the molecular hues trace a natal cloud core of GF 9-2. The highest contour of the H^^CO+ 
and N2H+ emission show east-west (EW) elongation. In addition, the N2H+ map shows a 
higher degree of central concentration than the H^^CO"'" and CCS maps with comparable 
beam sizes. The 3 mm dust continuum emission is located at the densest innermost region 
of the core, although its position is slightly displaced to the east of the H^^CO"*" and N2H+ 
peaks. 

Figure 5 shows velocity channel maps of the H^^CO+ (1^0) and N2H"'" (1-0) Fi,F=0,l- 
1,2 lines observed with the 45 m telescope; the latter transition is not blended with the other 
6 hyperfine components. The H^'^CO+ channel maps exhibit a velocity gradient along the 
major axis of the core; the redshifted gas dominantly lies to the NE of the 3 mm source, 
whereas the blueshifted one to the SW. 



3.2.3. Combined OVRO mm array and Nobeyama 45m Telescope Data 

Figure 6 shows the combined N2H"'' and H^'^CO^ images together with the original 45 m 
and OVRO maps, and Figure 7 shows channel maps of the combined data. We note that 
the OVRO and combined images are insensitive to the emission outside the fields of view of 
the OVRO array due to primary beam attenuation of the element antennas, and that the 
combined H^^CO"^ and N2H+ images are ~ 4 and 15 times less sensitive compared to the 
45 m maps, respectively. The OVRO images in Figure 6 trace the innermost region of the 
core and clearly show that both the N2H+ and H^^CO"^ peaks are shifted to the west of 
the 3 mm source. The combined images in Figure 6 show that there is a centrally peaked 
feature which is elongated along EW in the smaller scale and that intensity of the two lines 
decreases with increasing the radial distance from the center. However, contrary to Figure 5, 
one cannot see a well-defined velocity gradient in the combined channel maps of the H^^CO"'" 
emission (Figure 7). 
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3.3. Seeirch for Moleculeir Outflows and Jets 

3.3.1. Search for CO Outflow 

Figure 8 compares the blue- and redshifted ^^CO (1-0) and (3-2) emission integrated 
over the velocity ranges shown in Figure 2. In the figure, no highly-coUimated high- velocity 
"jet-hke" outfiow typically seen in class sources (Bachiller 1996, and references therein) 
is identified toward the 3 mm source, although formation of protostar in GF 9-2 core is 
strongly suggested by the H2O masers. This is because the terminal velocity of the blue- 
and redshifted emission, shifted by only ~ 3 km s^^ to Kjys (Figure 2), is too small to be 
interpreted as the line-of-sight components of CO outflow velocity (Vflow) in class sources 
(e.g., Vflow — 50 km s^^ for VLA 1623; Andre et al. 1990). Instead, a condensation of the 
low- velocity redshifted gas lies to the south of the 3 mm source in the (3-2) map. A similar 
condensation is marginally seen in the redshifted (1-0) map, but seems to be a part of the 
large-scale natal filament. These red condensations might represent a compact low- velocity 
outflow which has just launched from a protostar. On the other hand, a patchy distribution 
of the blueshifted (3-2) gas is observed, partly because of the quality of spectral baselines. 
The blueshifted (1-0) map shows a ridge to the NE of the 3 mm source, and is connected to 
the EW natal fllament, whereas the blueshifted (3-2) emission seems to be resolved into two 
peaks in the core center. 

First, we estimate an upper limit to the mass of the high-velocity outflow gas. Since 
the CO emission is very likely to be optically thick, we consider that the peak brightness 
temperature gives a good estimate of Tpx, assuming that the emission fllls the beam. We 
derive Tex of 5.9 K from the (1-0) and 6.9 K from (3-2) transitions, considering the cosmic 
background temperature (Tbg — 2.7 K). We attained 3(7 upper limits of 12 and 79 mK 
for the (1-0) and (3-2) transitions, respectively, with a velocity width of 50 km s~^ which 
is a typical velocity of high- velocity jet-like outflows associated with class sources (e.g., 
Bachiller 1996), although Figure 2 shows terminal velocity is ~ 3 km s~^. Given the mean 
value of T'ex= 6.4 K and [C0]/[H2] — 10'^ (Dickman 1978), we convert these upper limits 
into mass per lobe (Mobc) of 6.7 x 10"^ M© for (1-0) and 1.7 x lO"-'^ M© for (3-2). Clearly, 
these 3(7 detection limits in mass are sufficient to identify any of the CO outflows associated 
with class sources known to date (e.g.. Richer et al. 2000). 

Next, we estimate mass of the possible low-velocity (~ 3 km s~^) outflow suggested by 
the weak wing emission in Figure 2 and by the redshifted condensation in Figure 8. The 
extents of the lobes are approximately 0.053 and 0.032 pc for the blue- and red ones from 
the 3 mm source, respectively. Given Tex= 6.4 K, we obtain Miobe of 1.6 x 10^"^ Mq for the 
blue and 2.6 x lO"'^ Mq for the red. These lobe masses are considerably smaller than those 
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in class sources. 



3.3.2. Non-Detection of Shock Excited SiO Emission 

Thermal SiO emission is believed to be one of the most reliable shock tracers associated 
with molecular outflow (e.g., Codella, Bachiller & Reipurth 1999) because the destruction 
of dust grains releases silicon into the gas-phase and silicon is rapidly oxydized to produce 
SiO (e.g., Schilke et al. 1997). No evidence of outflow shocked regions is found in the SiO 
emission over a 0.22 pc x 0.14 pc region centered on the 3mm continuum source. The 
searched area is sufficiently large to detect such SiO emission enhanced at leading edge of 
the possible outflow lobe younger than dynamical age of <^ 0.053 pc/3 km s~^^ 1.7 x 10*^ 
years. This result suggests that the compact low-velocity outflow (§3.3.1) is not powerful 
enough to excite thermal SiO emission. 

4. Analysis 
4.1. Column Density of the Core 

4.1.1. N2H+ (1-0) 

To estimate the column density of N2H"'" (A^N2H+), we analyze the hyperfine structure 
(HFS) due to the electric quadrupole moment of the N nucleus. We performed least-squares 
fitting to each spectrum by applying a model spectrum having seven Gaussian components 
assumed to have the intrinsic intensity ratios in local thermodynamic equilibrium (LTE) at 
the limit of r = (Cazzoli et al. 1985). The details of the HFS analysis and the column 
density calculations are described in Appendix B.l. The free parameters are Tgx, velocity 
width (Ai'FWHM), LSR velocity of the brightest hyperfine (HP) component (J =1-0 Fi,F = 
2,3-1,2 transition at 93173.809 MHz), and the sum of the optical depths of the components 
(^tot)- We have presented an example N2H"'" spectrum toward the 3 mm peak position in 
Figure 3. We used only the single dish data because the velocity resolution is 3.4 times 
higher than that of the combined data. The parameters of the best-fit model profiles for the 
HF spectrum, shown by the green histogram, are summarized in Table 5. 

To perform HFS analysis on the combined data, we create a spectrum at each map pixel 
from the 3D cube data. The pixel size is set to be O'.'Q, the maximum size to avoid aliasing 
in the Fast Fourier Transformation (FFT) of the combined visibilities, whereas that for the 
45 m data is set to be 5". We obtained reasonable fitting results from 262 and 1776 spectra 
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for the 45 m and combined data, respectively. For the remaining spectra, we failed to achieve 
a fit because the weakest HF component of Fi,F = 0,1-1,2 whose intrinsic intensity is 0.111 
times of that of the brightest Fi,F = 2,3-1,2 did not have S/N-ratio higher than 5. 

Figure 9 shows total column density maps of N2H+ calculated from the 45 m and com- 
bined data. To produce Figure 9, we calculated A^N2H+ at each pixel using Eq.(B7) where 
we used Tgx, Ttot, and intrinsic velocity width, Avi^i, derived from the HFS analysis of the 
N2H+ spectrum at a pixel. The 45 m A^N2H+ map shows that the core has the densest region 
at the center, and that the column density decreases outward. The combined N2H+ map 
shows that the extent of the densest central region is smaller than the 45 m beam size. The 
region has an approximate extent of 10" (corresponding to 0.01 pc) and A^N2H+ of ^ 10^^ 
cm~^. Notice that the 3 mm source is shifted by 6" to the east of the A^N2H+ peak. 

To obtain fractional abundances of the probe molecules in §4.2, we recalculate iVN2H+ 
with various beam sizes. We smoothed the 45 m data by 40"and 80" beams for r < 20" 
and 40" regions centered on the 3 mm source. Table 6 summarizes the calculated A^N2H+ 
and compares CCS and NH3 column densities. We also reconstructed N2H+ images from the 
combined visibihties by giving the SHARC II beam as a synthesized beam size and calculated 
A^N2H+ inside the 3a level contour of the 350 /im emission. 

In addition, we analyzed the N2H"'" hue data with the 45 m telescope to use the hue as 
a temperature probe for the other molecules because the N2H+ HFS analysis gives rather 
accurate estimates of Tex and the spatial resolution for N2H+ is the highest among our data. 
We analyzed N2H"'" spectra inside the region enclosed by the half maximum, i.e., 50% level 

contour to the peak intensity and obtained mean Tex of 9.5±1.9 K (Table 5). We assume 
that the core has a temperature of 9.5 K and use the temperature to calculate the column 
density of the other molecules. Note that the adopted temperature is higher than the range 
of 5 to 7 K for the ambient gas inferred from the ^^CO (3-2) and (1-0) emission (§3.3.1). 



4.1.2. H^^CO+ (1-0) 

To estimate the column density of H^^CO^ (^Hi3co+)i we analyzed both the 45 m tele- 
scope and combined data. We firstly estimate optical depth (thi3co+) of the H^^CO^ (1^0) 
line toward the core center from an antenna temperature ratio of T^(H^^CO^)/T^(HC^^O"'"). 
Here, we assume that the beam filling factors, Tex, and the telescope beam efficiencies are 
the same for the two isotope lines. In addition, we assume that the solar abundance ratio 
pgjPg} of 5.5 is applied to the [Hi3C0+]/[HCi«0+] ratio. Using the peak T^ of 440 mK 
for the H^^CO"^ line in Table 5 (Tmb is converted into T^ by multiplying rfrah shown in Table 
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1) and the 1.5a upper limit of 39 mK in Tl for the HC^^O+ hne (Table 1), we found that 

the H^^CO"*" line is optically thin (< 0.1; Tabic 5). Given rHi3co+ of < 0.1 and T^x of 9.5 K 
(§4.1.1), we obtained iVHi3co+ from the 45 m and combined data (Table 6) by using Eqs.(B9) 
and (BIO) (see Appendix B.2), where the SHARCII beam was used as a synthesized beam 
to calculate A'^hi3co+ inside the 3a contour of the 350 /im emission. 

4.1.3. CsH2 2i2 — loi 

For calculating C3H2 column density (A^'c3H2) traced by the 2i2 — loi transition, we 
followed the calculations by Benson, Caselli & Myers (1998) (see their Appendix A2). We 
observed C3H2 only toward the core center. We estimated the optical depth using Eq.(17) by 
Benson et al. (1998) of Tc-^r.^ = In (^ y ^r'^-Tb )' ^CaHa = 0.34 for Tmb of 1.1 K, assuming 
Tex = 9.5 K (§4.1.1). We derived Av-mt (Table 5) from the second order moment of the 
spectrum because the line profile is not a single Gaussian (see Figure 2). Finally, a mean 
C3H2 column density is estimated to be (2.4 ± 0.3) x 10^"^ cm^^ (Table 6). 

4.1.4. CCS A3 -32 

In the similar manner as the H^^GO^ analysis, we estimated the GGS column density 
(iVccs) using Eq.(Bll). We obtain rccs of 0.94 (Table 5) from aratio of TX(GGS)/TX(GG34S) 
= 15.2 which is derived from the peak of 560 mK for GGS (Table 5) and a 1.5a upper 
limit of 37.5 mK for GG'^'^S (Table 1). Here we adopt the solar abundance ratio p^S]/p^S] of 
23. Given rccs of 0.94 and Tex of 9.5 K, we obtain Nccs of (1.5 ± 0.2) x 10^^ cm-^ toward 
the core center with a 40" beam (Table 6) . 

4.1.5. NHs (1,1) 

We present spectral line profiles of the NH3 (J, K) = (1, 1) and (2,2) lines in Figure 10. 
The (1,1) spectrum consists of five groups of the lines. The most intense emission in the 
(1,1) transition, the central group which contains eight of HF fines, has a peak of 640 
mK, whereas no (2,2) emission is seen above our 3a level detection threshold (Tmb = 210 
mK). Neither the (3,3) nor (4,4) emission was detected. 

To estimate NH3 column density (A'^nhs), we followed the previous works, e.g., Win- 
newisser, Ghurchwell & Walmsley (1979), Ungerechts, Walmsley & Winnewisser (1980), 
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Stutzki, Ungerechts & Winnewisser (1982), and Ho & Townes (1983). We analyzed the 
spectra taken inside the 50% level contour of the total intensity map (Figure 4d). In the 
similar manner as the N2H"'' HFS analysis, our ammonia HPS analysis gives Ttot and Avfwhm, 
and the LSR- velocity of the strongest HF component. The brightest HP emission has optical 
depth of ~ 1.6 toward the core center (Table 5). Here wc assume that the beam filling factor 
is unity. It is known that rotational temperature (Tr^2i) between the (1,1) and (2,2) tran- 
sitions and the kinetic temperature (T^) agree well with each other for high density 10^ 
cm~^) and cold (^ 10 K) gas (Tafalla et al. 2004). We therefore assume Tj.^2i = ^k, and 
that Tk is equal to Tex(N2H+) of 9.5 K (§4.1.1). Now we calculate a mean column density of 
(1.6 ± 0.6) X 10-*^^ cm"^ for the NH3 molecule in the (1,1) state inside the 50% level contour 
of the total integrated intensity map by inverting Eq.(2) of Ungerechts ct al. (1980) Wc 
consequently obtained total NH3 column density from the (1,1) column density by Eq.(2) of 
Winnewisser et al. (1979) (Table 6). 



4.2. Fractional Abundances of the Molecules in the Core 

To calculate mass (Mlte) of the core from the column densities (§4.1), we need to find 
the fractional abundances (X^) of the probe molecules to H2 molecules. Such an abundance 
is known to vary under different physical and chemical conditions. Bergin & Langer (1997) 
theoretically predicted that all molecules are depleted in some degree from gas phase for high 
densities of n(H2) ^ 10^ cm~^, regardless of dynamical evolutionary phases. The depletion 
of sulfur-bearing molecules such as CS, SO, and CCS is particularly sensitive to the density 
increase because they tend to stick to both H2O and CO ice covered grains. In contrast, 
nitrogen-bearing molecules of NH3 and N2H+ do not seem to be susceptible to freezing out 
on the grain surface (Bergin & Langer 1997). This prediction is confirmed by the recent 
observations toward starless cores, e.g., Tafalla et al. (2002) who showed that N2H"'" and 
NH3 column densities correlate better with the dust column density than molecules such as 
CO and CS. This characteristic, coupled with high angular resolution, makes N2H+ ideal 
as a reference molecule. In addition, we believe that the possible low-velocity compact 
outflow does not affect our discussion about molecular abundances, considering our flndings 
from radial density profile, position-velocity diagrams and velocity width maps of the core 
(described in §4.4 and §4.5). 
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4.2.1. N2H+ and H^^CCP- 

The abundances of N2H+ and H^'^CO^ arc estimated through comparison to the mean H2 
column density, A^(H2)g//, derived from the 350 fim continuum flux with the SHARC II beam 
(§3.1.2). Because the extent of the 350 /xm emission is approximately 4 times smaller than 
those of the N2H+ and H^^CO+ emitting regions, such a comparison is possible only in the 
common innermost region of the core. Comparing the A'"(H2)g// with the N2H''" and H^^CO"^ 
column densities with the SHARC II beam size (Table 6), we obtained X(N2H+)g// =(9.5 ± 5.0) x 10" 
and X(Hi=^CO+)g// =(8.5 ±4.7) x 10"^^ (^able 7). Note that the errors in X^n are mostly 
due to the uncertainties in (§3.1.2). 

We argue that N2H+ in GF9-2 is not depleted, whereas H^^CO+ is depleted. Our 
estimate of X(N2H+)g// is comparable to those in the following studies. Womack et al. 

(1992) obtained a typical value of X(N2H+)= 7 x 10"^° in the dark cloud cores TMC-1 and 
L134N (^HPBW = 70"). Benson et al. (1998) derived (7 ± 4) x 10"^° from a Haystack 37m 
(22") survey towards 60 low-mass cores. Casclli et al. (2002) reported (3 ±2) x 10^^° from a 
FCRAO 13.6 m (54") mapping survey towards 34 dense cores. Tafalla et al. (2002) reported 
the range of (0.75 - 2.4) X 10^1° toward 5 starless cores using the FCRAO. The consistency 
between our measurements with 8" beam and those with 22", 54", and 70" beams suggest 
that N2H+ in GF9-2 is not depleted at the core center. On the other hand, X(H^^CO+)g// in 
GF 9-2 is by a factor of almost 6 smaller than the lowest abundances in other similar cores 
measured with 67" beam (Table 9 by Butner et al. 1995), e.g., 5.6 x 10"^^ for B335 (class 0) 
and 5.4 x 10"" for L 1495 (class I). This fact suggests that H13qo+ in GF9-2 is very likely 
to be depleted at the center, although we need to verify the radial variation of X(H^^CO"'"). 



4.2.2. C3H2, CCS and NH3 

Using the X(N2H+)g// as a reference value, we estimated the abundances of C3H2, CCS 
and NH3. For this purpose, we assumed that N2H"'" in GF 9-2 docs not have radial abundance 
variation. Recall that Tafalla et al. (2002) reported that N2H+ abundances are spatially con- 
stant in the 5 starless cores they surveyed. Given the X(N2H"'")g// and a column density ratio 
of N(C3H2)/N(N2H+) with 20" beam (Table 6), we obtained X(C3H2)2o'/ of (2 ± 1) x lO'^o 
(Table 7). Similarly we obtained X(CCS)^q// of (5 ± 3) x 10-^'' from an N(CCS)/N(N2H+) 
ratio for r < 20", and X(NH3)gQ// of (2 ± 1) x 10"^ from an N(NH3)/N(N2H+) ratio for 
r < 40". 

The X(C3H2)2q" value in GF9-2 would fall at the lower end of the X(C3H2) range, 
(3.2-40) X 10"^^ obtained from the survey toward dense dark clouds with 54" beam (Cox, 
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Walmsley & Giisten 1989), although the beam sizes differ by a factor of ~ 2.5. The estimated 
X(CCS)^p// = (5 ± 3) X 10""^° is consistent with those of the starless cores measured with 
comparable beam sizes, i.e., 5.8 x 10"^° in TMC-ID and 4.9 x 10-^° in L 1498 with 45" -53" 
beams (Wolkovitch et al. 1997), although smaller than (8 ± 3) x 10~^ with ~ 52" beam in 
L1521F (Shinnaga et al. 2004). On the other hand, the X(CCS)^q// in GF9-2 is found to 
be larger than those in the representative class sources B 1 and B 335 (~ 9 x 10^^^ and 
~ 1.6 X 10"^°, respectively) reported in Lai & Crutcher (2000) with ~ 30" BIMA beam. 
Lastly, the X(NH3)gQ// = (2 ± 1) x 10""^ in GF9-2 is somewhat smaller than the mean value 
derived from the large compilation of literature toward the cores in the Taurus and Ophiucus 
regions 

((4.0lo:8) X 10"^ Table BIO of Jijina et al. 1999), but is comparable to the range of 
(0.4 - 2.3) X 10"^ for the 5 starless cores surveyed by Tafalla et al. (2002). 

Here we wish to discuss age of the GF 9-2 core in terms of chemical evolution. Suzuki 
et al. (1992) analyzed their 45 m telescope CCS survey data by developing a chemical evolu- 
tionary model, and concluded that CCS is the most abundant in the early evolutionary stage 
for a quiescent dark cloud core, whereas NH3 is less abundant and becomes abundant in the 
later evolutionary stage. Aikawa et al. (2001) advanced such a model for a dynamically 
collapsing core and pointed out that N2H"'" shows similar chemical evolution as NH3. Com- 
paring the X(NH3)gQ// and X(CCS)^q// with Figure 16 in Suzuki et al. (1992) in conjunction 
with Figure 5 of Aikawa et al. (2001), we argue that GF9-2 core has experienced 10^ — 10^ 
years since the diffuse stage. Because of the large uncertainty in our abundance estimate, 
we have limited our discussion to the 3 molecules whose abundance variation with time are 
large enough to compare with our results. 



4.3. LTE Mass of the Core 

Using the column densities of the N2H+, R^^CO+, CCS, and NH3 lines (§4.1) and their 
fractional abundances (§4.2), we now calculate the LTE mass (Mlte) of the gas in the area 
(A) enclosed by the 50% level contour of each total integrated intensity map of Figure 4 (see 
Table 8). The mass is given by Mlte = A*?71h • ' ^ where /i is the mean molecular weight 
(2.33 for [He]= 0.1[H]) and mn is the atomic mass of hydrogen. The obtained LTE masses 
range from 1.2 to 4.5 Mq. 
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4.4. Radial Column Density Profile of the Core 

We derive radial column density profile of GF 9-2 core using the N2H+ column density 
maps. Recall that the N2H+ in GF 9-2 is unlikely to be depleted over the full extent of the 
core, whereas the H^^CO+ may be depleted (§4.2.1). We assume that the core is spherical. 
To azimuthally average the column density maps, we take the N2H+ emission peak as the 
center. Figure 11 represents the resultant column density profiles. The profiles can be well 
fitted by a power-law function of N(r) — A^o(^)~^ where p is the power-law index, and A^o 
is the column density at r = ro = 1". By minimizing the values, we obtained the best- fit 
parameters of j9 = 1.0 ± 0.1 and A^o = (7.3 ± 1.1) x 10^^ cm~^ for the combined data, and 
p = 1.1 ± 0.2 and Nq = (8.5 ± 2.5) x 10"*^^ cm~^ for the 45 m telescope data. The two fitting 
results show good consistency within uncertainties, indicating that a negligible amount of 
the dense gas has been evacuated by the possible compact outflow in the innermost region 
of r ^ 3". 



4.5. Velocity Structure of the Core 

4.5.1. Overall Velocity Structure 

To examine the velocity structure of GF 9-2 core, we made position- velocity (PV) di- 
agrams using the 45m N2H+ (1-0) Fi,F=0,l-l,2 and Hi3C0+ (1-0) line data (Figure 12). 
We also overlaid PV diagrams made from the combined H^^CO+ data set. Since the core 
is elongated along the NE-SW (Figures 4, 5 and 9), we select the two cutting axes at P.A. 
of 45°and 135°, passing through the position of the 3mm source as the major and minor 
axes, respectively. Along the major axis, the 45 m data show a weak velocity gradient from 
NE toward SW in both the molecular lines. The gradient along the major axis is an order 
of 1 km s~^ pc~^, comparable to the typical value measured in the Taurus dark cloud cores 
(Goodman et al. 1993). If the gradient represents the rigid-body-like rotation of the core, it 
is converted to an angular velocity of A ~ 3 x 10~^^ s~^. On the other hand, the minor axis 
plots do not show a well-defined velocity gradient commonly seen in the two lines, although 
the H^^GG"*" diagram appears to show an SE-NW gradient if we concentrate on the lowest 
two contours. 

Another important finding from Figure 12 is the presence of broader line width region 
at the core center, prominently seen for the H^^GO"*". The blueshifted gas is seen up to ~ 0.6 
km s~^ from the Vgys at the central region of r ^ 5" (corresponding to 1000 AU), whereas 
the redshifted counterpart about ^ 0.3 km s"-^. The presence of such blueshifted emission 
and the weak redshifted counterpart at the core center (see also Figure 2) suggests that gas 
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is infalling where redshifted gas is obscured by the foreground cold gas along the line-of-sight 
(e.g., Walker et al. 1986). Notice that the line broadening at the core center is unlikely to 
be caused by the possible compact outflow because the red lobe is more intense than the 
blue one (see Figure 8), inconsistent with the above results. 

4.5.2. Velocity Width 

We firstly estimate intrinsic velocity width {Avi^t) toward the core center by analyzing 
the N2H+, C3H2, CCS, and NH3 fines (Table 5). Note that Avi^t is the velocity 

width corrected for line broadening due to instrumental velocity resolution and optical depth 
of the line (Phillips et al. 1979). Both corrections are made for the C3H2 and CCS lines. 
The latter correction is not made for N2H+ and NH3 data because HFS analysis obtains the 
optical-depth-free AfpwHM (see Appendix B), and not for the H^^CO+ data because it is 
optically thin. The obtained Avint toward the center is between 0.35 and 0.61 km s~^. This 
range is 2~3.5 times larger than the isothermal sound speed (cg) of the ambient gas of 0.18 
km s-i at Tk = 9.5 K (§4.1.1). 

Subsequently, we analyzed velocity width structure of the core. Figure 13 shows the 
Avint maps made by spectral moment analysis of the 45 m telescope N2H+ (1-0) Fi,F=0,l- 
1,2 and H^^CO+ (1-0) fine data. We adopt the moment method to keep a consistency in 
the analysis between the two lines. To make the maps, we made spectra from the cube data 
at each point, then searched for a peak LSR- velocity in each spectrum to determine two 
velocity ranges from the peak velocity toward both the blue- and redshifted sides where the 
intensities are above the 1.5a level. We calculated the second-order spectral moments over 
the velocity ranges defined for the individual spectra. For the N2H+ map (Figure 13b), we 
corrected the line broadening effect due to optical depth. We used Eq.(3) in Phillips et al. 
(1979) by substituting Tp of Eq.(3) with rtot x s where s is the intrinsic relative intensity of 
the Fi,F = 0,1-1,2 transition (0.111; Cazzoli et al. 1986). The maximum line broadening at 
the most opaque (rtot x 0.111 = 1.5) position is 1.3 times the width for Tp — > 0. Considering 
the thermal noise and the velocity resolutions, we estimate uncertainties of 0.11 and 0.12 
km s~^ in Av^^^t for N2H+ and H^^CO"*", respectively. As seen in PV diagrams (§4.5.1), both 
the N2H+ and H^^CO^ maps show a clear increase of Afint toward the core center. The 
peak of the broadening does not fall onto the core center. It also should be noted that Av[^i 
of the core gas exceeds Cg of 0.18 km s^^ over the core (sec also Table 5). Furthermore, 
comparing Figures 8 and 13, one cannot discern any clear effects such as an evacuated shell 
by the possible compact outflow in Figure 8; the ^^CO (3-2) condensations are only observed 
toward the 3 mm continuum source. 
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The above analysis clearly shows that there exist non-thermal motions in the core. 
Wc estimate non-thermal velocity width (Avnth) by subtracting the thermal velocity width 
(Avthm) from Ai>int: (At>int)^ = (A^thm)^ + (Avnth)^- Here the thermal contribution is 
given by Avthm = . /5i!£-?MZk where k is the Boltzmann constant, and mmoi is the molecular 

y '^mol 

weight: 29 amu for N2H+, 30 amu for H^'^CO^, 56 amu for CCS, and 17 amu for NH3. 
Assuming that is equal to T^x of 9.5 K (§4.1.1) and that is constant over the entire 
core, At;thm is calculated to be 0.12 km s"^ for N2H+, 0.12 km s'^ for Hi3C0+, 0.09 km s"^ 
for CCS, and 0.16 km s~^ for NH3 lines. This leads to mean non-thermal velocity widths 
over the core of (Avnth) = 0.21, 0.58, 0.40, and 0.46 km s"^ for the N2H+, E.^^CO+, CCS, 
and NH3 lines, respectively, suggesting the presence of supersonic gas motions. In obtaining 
(Ai>nth), we used the (Ai>int) values in Table 5. 



5. Discussion 

5.1. Evolutionary Stage of the GF 9-2 Core: A Very Young Protostar Before 

Formation of a Lcirge-Scale Moleculcir Outflow 

With its low luminosity (~ 0.3 Lq) and cold SED (Tboi < 20 K; Wiesemeyer et al. 
1997), we believe that the protostar embedded in GF9-2 core is at an extremely early 
phase of low-mass star formation, evidenced by the absence of a large-scale high-velocity 
molecular outflow, the possible detection of a compact low-velocity outflow, and the detection 
of H2O masers. We estimate the protostar's age from dynamical age (tdyn) of the possible 
compact low- velocity outflow (§3.3.1). The terminal velocity (Vt) of the spectra would give 
the line-of-sight component of the outflow velocity of |T4 — V^ysl — ^owcosi where i is the 
inclination angle of the outflow axis to the line-of-sight. Given the projected lobe length of 
ZiobeSini = 1, 1000 and 6600 AU and Vflowcosi ~ 3.3 and 3.1 km (§3.3.1), we have idyn = 
hohe/Vaow ^ 1-6 X 10^ and 1.0 x 10^ years with i = 45° for the blue- and rcdshiftcd ^^CO (3-2) 
lobes, respectively. These dynamical ages are comparable to those for the extremely young 
class sources, e.g., (0.6-3) x 10^ years for VIA 1623 and (0.8-2) x 10^ years for IRAM 04191 
(Table 1 in Andre, Motte & Backman 1999, and references therein). 

The youth of GF9-2 is demonstrated in Figure 14 where we present a plot of CO 
outflow momentum rate (Fcq) normalized by the bolometric luminosity (-Fco ' c/Lboi) vs. 
Tboi for nearby (d < 350 pc) low-mass (Lboi < 50 Lq) YSOs. We employ Fco-c/Lboi 
because molecular outflows in low- mass protostars are likely to be momentum-driven (e.g., 
Cabrit & Bertout 1992; Richer et al. 2000, and references therein). The temperature Tboi 
is believed to trace the evolutionary stage of YSOs (Myers & Ladd 1993). We calculate 3a 
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upper limit of Fco • c/Lboi = ■ 7^ ~ 9 and 20 for the ^^CO (3-2) blue- and red 

Hobe ^bol 

lobes, respectively, with i — 45°. Figure 14 illustrates that "class 0" is the most powerful 
outflow phase and "class I" is the decline phase (e.g., Moriarty-Schicvcn et al. 1994). The 
most remarkable result is that GF 9-2 shows the lowest outflow activity, among the coldest 
low-mass class sources, although there exists enough circumstellar material (§4.3) to be 
dragged by a possible protostellar jet to develop a large-scale molecular outflow. Note that 
the prototypical class sources of VIA 1623 (AWB93) and IRAM 04191 (Andre et al. 1999) 
have powerful CO outflows, despite that their bolometric luminosities are comparable to that 
of GF9-2: 1 for VLA1623 and 0.15 Lq for IRAM 04191. We thus believe that the low 
luminosity of GF 9-2 is not responsible for the lowest outflow activity. Gonsequently, it is 
very likely that the protostar of our interest is the youngest class protostar ever identified 
and is just at the beginning of the outflow launch. 

5.2. The Initial Conditions for Gravitational Collapse of a Low-Mass Star 

Forming Core 

5.2.1. Stability of the GF9-2 Core 

If we interpret that the non-thermal velocity width obtained in GF9-2 core (§4.5.2) is 
attributed to turbulent motions, we expect that the turbulence can help support the core 
against its self-gravity. To test this hypothesis, we compare the LTE mass of the core with 
virial mass Myir including the turbulence. Assuming a spherical system, the virial mass 
is given by Myir ~ it!eff(Ai'int)^/G' where (Ai>int) includes the non-thermal width, i.e., the 
turbulence. Given (Ai>int) and R^s values (Tables 5 and 8), we found that M^-„ ranges from 
0.39 to 2.9 M0(casc 1 in Table 8). Since Mvir — Mlte within the uncertainties, the core is 
most likely to be in virial equilibrium. However, it should be noted that this hypothesis does 
not reconcile with the line broadening toward the core center (Figure 13). Moreover, as far 
as we know, no cores which are supported by turbulence exhibit a p(r) oc r^^ density proflles 
(see Figure 11). Therefore, we conclude that the non-thermal line width is not caused by 
turbulence. 

Alternatively, if we consider that infall causes the non-thermal line width, that is, if 
GF 9-2 core is gravitationally unstable, the drawbacks in the previous hypothesis would be 
solved. The presence of infalUng motion is supported by the following facts: (1) the presence 
of the blueshifted gas and the less prominent redshifted counterpart in the central r ^ 1000 
AU (§4.5.1), (2) the spectroscopic signature of infall, i.e., the "blue-skewed" proflle observed 
in the optically thick N2H+ (1-0) Fi,F=0,l-l,2 line (Figure 2), and (3) the detection of 
the H2O masers, a clear signpost of star formation, naturally leading to the presence of 
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accretion flow onto the star-disk system. We, however, point out a drawback in the above 
(2); the optically thin H^^CO"'" emission docs not show a single Gaussian profile expected 
from the infall interpretation, while the CCS does. In this interpretation, Myir should be 
estimated only from the thermal motion as M^i,. = The estimated M^i^. from the 

4 lines takes a range between 0.078 and 0.25 Mq (case II in Table 8) which is one order of 
magnitude smaller than Mlte- The small Myir is consistent with our interpretation that the 
core is under gravitational collapse. As discussed in §4.5.2, the core would be rotating with 
A ~ 10~^^ s~^. The collapse cannot be supposed by the centrifugal force due to the rotation 
because the rotational energy is less than a few percent of total energy required to support 
the core against its self-gravity. For such an unstable core, the gas is expected to be more 
perturbed toward the center, and the velocity width should increase inward as a result. In 
fact. Figure 13 clearly indicates that Avint in GF9-2 increases up to ~ 0.6 km s~^ toward 
the center. 



5.2.2. Initial Conditions for the Gravitational Collapse of the Core 

To bring out the initial conditions for the gravitational collapse in GF9-2 core, we 

examine the density and velocity structures by comparing with some key characteristics of the 
Shu (1977) and Larson-Pcnston (Larson 1969; Pcnston 1969) solutions. The radial H2 column 

density profile of the core is given by (7.7±4.1) x 10^3 (^^^iji^^"^ cm'^ for 

0.001 < ^ < 0.05 from the combined data and (8.9±4.5) x lO^^ (|g^2S^ j ^ (2^X11) 
cm~^ for 0.005 < ^ < 0.08 from the 45 m data (Figure 11). It should be noted that the dense 
N2H+ gas seem to obey the power-law profile with the index of —1 even at the innermost 
region of r ^ 1000 AU where a possible compact outflow might begin to destroy the core. 
The power-law indices of —1.0 and —1.1, converted into —2.0 and —2.1 for the volume 
density, well agree with those in the two solutions after the formation of the protostar (Shu 
1977; Whitworth & Summers 1985). The best-fit H2 column densities at r = 200 AU fall 
between the two extreme predictions at T = 9.5 K: (2.8 ± 1.5) x 10^^ cm~^ for the Shu and 
(13 ± 7) X 10^3 cm~^ for the Larson-Penston solutions. Here we used a relation of iVo(H2) = 
1.0 and 4.4 X 2^ ■ ^^^^^^ for the Shu's and Larson-Penston's cases, respectively (Hunter 
1977; Shu 1977; Foster & Chevalier 1993). Because of the abundance uncertainty, however, 
we cannot conclude which solution describes well the collapse in GF 9-2 core. 

Contrary to the above discussion based on the density profile, the velocity width maps 
(Figure 13) help to discriminate the two theoretical scenarios for the collapse of GF 9-2 core. 
Notice that the infiuence of the compact outfiow on the velocity width cannot be discerned 
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toward the 3 mm source. Figure 13 shows that the velocity width decreases outward, but 
remains ~ 0.3 km s~^, which is twice the ambient Cg, even in the outermost region of r ~ 50" 
(corresponding to 0.05 pc) where the density profile shows p(r) oc (Figure 11). Recall 
that p(r) oc profile characterizes outer part of a collapsing core. These facts clearly 
indicate that the extended Larson-Penston solution for t > is more likely to be the case 
because it requires the infalling state when a protostar is about to be born (f (r) ~ — 3.28cs 
at i — > —0; Hunter 1977), whereas the Shu solution requires a static state {v{r) — a,t t — 
0). It is, therefore, highly probable that the gravitational collapse of GF9-2 core is well 
described by the runaway collapse scenario. 

Subsequently, we estimate the age of the GF9-2 protostar (tprotostar) based on the lack 
of the density profile of p{r) oc r~^/^ down to the innermost region of r ~ 3" (see Figure 11). 
Recall that p(r) oc r~^/^ profile is expected for a central free-fall region which appears in both 
of the solutions at i > 0. The non-detection of a free-fall density profile indicates that the 
central free-fall region, propagating outward with supersonic velocity, has not yet developed 
up to a radius (rg) of ~ 600 AU, suggestive of the youth of the GF9-2 protostar. However, 
we do not exclude the possibility that the absence of a detection of the p(r) oc r^'^^^ profile 
at the center could be caused by the formation of binary (described in §5.3) whose projected 
separation is ~ 1200 AU. Using Eqs. (9) and (10) in Ogino, Tomisaka & Nakamura (1999), 

we estimate tprotostar to be less than tprotostar — 5000 ( g ^ years since the formation of 
the central protostar for the Larson-Penston solution. Although the age estimate is a factor 
of 2-3 smaller than the outflow dynamical age of tdya — (1-0 — 1.6) x 10^ years (§5.1), we 
argue that the central protostar has just begun to exhibit its jet activity traced by the H2O 
masers (§5.1). 

Consequently, we propose that GF9-2 core has started its gravitational collapse with 
initially unstable state. Here, the term initial means t — —tg where ts is free-fall time of the 
core (Larson 1969; Penston 1969). For GF9-2, we calculate iff ~ 2 x 10^ years. We therefore 
believe that the age of the core, tcore = ^ff + ^protostar, would be ~ 2 x 10^ years, which has 
a consistency with the time scale of chemical evolution (§4.2.2). 



5.2.3. Infall Velocity and Mass Accretion Rate in the Core 

Assuming that Avmt represents a velocity difference of blue- and redshifted infalling gas 
along the hne-of-sight, let us discuss infall velocity (vint) of the core gas. We found that Vinf 
varies from 0.2 to 0.4 km s"^ using a relation of vinf ^ llfinf.red — 'J^inf,biue|/ cosi ~ 
the Avint range of (0.28 - 0.63) km s"-*^ over the core (see Figure 15) and i = 45°. We now 

calculate the mass infall rate from ^ = 47rit!gg//mH •?T'H2(-Refr) •'^inf — '^T^ReS • ^ "x(N2H+) ""^^ ' 
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Vinf for a sphere with the radius of i?eff- Here nH2(-Refr) is a volume density of a thin shell 
whose radius is Res and is approximated as 2j^ff .x(N2H+) " ^i"^®^ '^mf — 0-29 km s"^ at r ~ Res 
from the N2H+ Av-mt map (Figure 13), we have ^ ~ 3 x 10~^ Mq year~^ with uncertainty 
of a factor 1.6 for i?cff in Table 8, the column density profile (§4.4), and in Tabic 7. The 
inferred mass accretion rate falls within the range between the predicted values of 1.5 x 10^® 
and 7.4 x 10^^ Mq ycar^-*^ for the Shu and Larson-Penston solutions, respectively, at 10 K. 
These accretion rates arc calculated by mocf^^/G where itlq is the non-dimensional factor of 
0.975 for the Shu and 46.92 for the Larson-Penston solutions (Hunter 1977). Our small mass 
accretion rate, compared to that of the Larson-Penston's solution, is attributed to both the 
smaller column density and infall velocity. 



5.3. Is a ProtobinEiry System Forming ? 

Wc found that the position of 3 mm continuum source is displaced by ~ 6" to the cast of 
the N2H+ column density peak, which would be explained by the presence of a protobinary 
system. In this interpretation, the 3 mm source represents a protostar and the N2H+ peak 
a preprotostellar condensation. This is because the N2H+ condensation is deeply embedded 
in the extended 350 //m emission (Figure 9b), but is not accompanied with a compact 3 mm 
source. To investigate the nature of the binary members, we made a Avint map from the 
combined H^'^CO^ data shown in Figure 15 as the line is expected to be optically thin. In 
§3.1.2, wc asserted that the 3 mm contimmm emission would represent circumstellar material 
around a protostar. The line broadening toward the 3 mm source seen in Figures 12 and 15 
suggests the local enhancement of mass accretion, leading to the existence of a protostar. 
On the other hand, the N2H+ peak (white contour) seems to be a gas condensation or a 
preprotostellar core rather than an apparent peak caused by local abundance enhancement 
because X(N2H+)g// inferred over the 3a level contour of 350 //m emission is comparable 
to those obtained with low spatial-resolution studies (§4.2.1). It should be noted that the 
H^^CO"*" emission has a secondary peak near the N2H"'" condensation as seen in the OVRO 
maps (Figure 6). Furthermore, the combined N2H"'" map (Figure 6) does not show a shell- 
like structure which is expected to be seen if N2H"'" disappears from the gas phase at the 
core center (e.g., Kuipcr, Vclusamy & Langer 1996; Belloche & Andre 2004). Wc therefore 
conclude that the positional displacement of N2H+ peak is not caused by chemical evolution, 
but the N2H"'" peak represents a preprotostellar condensation, making it younger than the 
eastern object. The detection of the two objects reinforces the assertion that we witness a 
forming binary system with a separation of ~ 1200 AU. 

The masses of the circumstellar material (Mem) for each binary member would be an 
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order of 0.1 Mq. For the eastern object, the 3 mm continuum flux density (§3.1.2) integrated 
inside its 3a level contour (see Figure 15) implies that the gas plus dust mass ranges from 
~ 0.1 (case I in §3.1.2) to 0.2 Mq (case II). For the western object, we estimate (0.1 - 0.3) 
Mq from the N2II'*" intensity integrated inside the white contour, whose extent is similar 
to that of the 3 mm source (Figure 15). The inferred Mem seem reasonable for an amount 
of circumstellar material with 100 AU scale around a newly formed protostar at the very 
early accretion phase because of Mem ^ protostar = ^ ^ ^protostar — 0.15 Mq (§5.2.2 and 
5.2.3). In this context, the EW elongated structure at the core center seen in the 45 m maps 
(Figure 6) may be interpreted as a 1000 AU scale circumbinary disk-hke structure because 
the size would be smaller than the extent of the protostellar envelope traced by the 350 /im 
continuum emission (6600x4500 AU; §3.1.2). It is impossible to uniquely define spatial 
extent of the disk-like structure. However, if we arbitrarily define the disk-like structure by 
a region enclosed by the 18 or 21a level contours of the combined H^^CO"'" map (Figure 15), 
mass of the circumbinary disk-like structure (Mcbd) calculated from the N2H"'" data would 
range (0.2 - 0.5) Mq with an approximate effective radius of ~ 2000 AU. Considering the 
uncertainties in spatial extents and in kq, P (§3.1.2) and X(N2H+) (§4.2.1), the mass of disk- 
hke structure is likely to be plausible to be hosted in the (0.4 - 0.8) Mq envelope (§3.1.2) 
and to host the two binary members with the 0.1 Mq order circumstellar material. The 
mass distribution of Mcorc > Meny ^ -^cbd ^ Afcm does not contradict with our hypothesis 
that the central protostar (s) is at extremely early evolutionary phase. 

If a circumbinary disk exists, the radius of disk-like structure would be physically related 
to the rotation of the natal core, described in §4.5.1. Although the axis for the large-scale 
rotation of the core (P.A.~ 135°) is not perpendicular to the major axis of the disk-like 
structure (P.A.^ 90°), a rough estimate of radius of the disk-like structure would be given 
by the centrifugal radius, /?ccnt) calculated from the specific angular moment of the core 
gas. With a relation of 7?ccnt = (4i?efr^^)/(<^^cbd), we obtained R,.„,t - 900 AU for A of 
3 X 10~^^ s~^ (§4.5.1) and i?efr in Table 8. The centrifugal radius is roughly twice the radius 
of disk-like structure, probably supporting our interpretation. Lastly, let us estimate the 
Toomre's Q parameter to examine whether or not gravitational fragmentation of the disk-like 
structure led to the formation of the binary. The Q parameter is given by for Keplerian 
rotation. Here S is the surface density of the disk-like structure given by Mcbd/-4.cbd; an 
area of ^cbd — 4 x 10^^ cm^ is calculated for the region described above. We obtained 

Q - 0-01 (s^om-r) (wk)'" {^) ( dxi^'cmO - Considering the uncertainties in the 
input parameters, it is certain that the Q value is less than unity. We, therefore, suggest 
that the binary system is likely to be forming in the core center due to the gravitational 
fragmentation of the circumbinary disk-like structure, whose radius appears to be determined 
by the rotation of the core. 
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6. Summciry and Future Works 

The main findings of tliis work and future considerations are summarized as follows. 

1. Our high density gas observations identified a molecular cloud core with a size of ~ 
0.08 pc toward GF9-2. In the core center, a protostellar envelope with a size of ~ 
4500 AU is identified through 350 /xm continuum emission. The mass of the core and 
the envelope are ~ 3 and ~ 0.6 Mq, respectively. Obtaining fractional abundances 
of the probe molecules, we propose that GF 9-2 core should be at an early phase of 
chemical evolution, although the radial variations of the abundances should be revealed 
to accurately assess the chemical evolutionary stage. 

2. Our deep search for CO wing emission revealed that the protostar in GF9-2 has not 
developed a large-scale high-velocity "jet-like" outflow. Instead, a compact low- velocity 
outflow may be formed at the core center. Thus, GF9-2 is very likely to be at an 
extremely early stage in low-mass star formation before arriving its most active outflow 
phase. To test the conclusion, we need to conduct high angular resolution and deeper 
search for compact outflow(s) in higher- J transitions of CO. 

3. The radial column density profile of the core can be well fitted by a power-law form of 
p{r) oc for 0.003 < r/pc < 0.08 region. The power-law index of —2 agrees with the 
expectation for an outer part of gravitationally collapsing core. On the other hand, 
based on the absence of a p(r) oc r~^^^ profile, we found that a central free-fall region 
has not developed in a region larger than r ~ 600 AU. 

4. The PV analysis identified a large-scale velocity gradient with an order of 1 km s~^ 
pc~^ along the major axis of the core. The gradient is interpreted as a rigid-body-hke 
rotation with an angular velocity of ~ 3 x 10~^^ s~^, comparable to those measured in 
other similar cores. We found that the gas is supersonic all over the core and that the 
velocity width is enhanced up to ~ 0.6 km s^^ in the central 1000 AU toward the 3 mm 
source. The supersonic velocity width and the p(r) oc r^^ profile are explained by infall 
in the extension of the runaway collapse scenario proposed by Larson and Penston for 
t > 0. The current mass accretion rate of ~ 3 x 10~^ Mq year~^ is estimated from 
^;inf ~ 0.3 km s~^ at r ~ 7000 AU. We discuss the likehhood that GF 9-2 core has 
been undergoing gravitational collapse ^ 5000 years, starting from initially unstable 
conditions. 

5. We discovered a potential protobinary system with a projected separation of ~ 1200 
AU embedded in a circumbinary disk-like structure with ~ 2000 AU radius at the core 
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center. The binary consists of a very young protostar and a preprotostellar condensa- 
tion. The circumbinary material mass is on the order of 0.1 Mq. Taking account of the 
centrifugal radius calculated from the core rotation and the Toomre's Q value for the 
disk, we suggest that the binary formation may have been initiated by gravitational 
fragmentation of the circumbinary disk-like structure. Further studies are needed using 
high-angular-resolution deep imaging in submm and FIR continuum emission. 
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A. Combining the OVRO mm-array and Nobeyama 45 m Telescope Data 

We describe a method of combining the OVRO mm-array and Nobeyama 45 m telescope 
data to obtain high-spatial-resolution and missing-fiux free images. We used the 45 m data to 
fill the central u-v hole of the mm-array visibility data. Since the OVRO array consists of six 
10.4 m diameter telescopes, we can attain the minimum projected baseline length normalized 
by the observing wavelength (-D™|°t) of ~ 3.5 k\ at 3 mm when the array takes the most 
compact "C" configuration. The minimum spatial frequency is sufficiently smaller than the 
maximum spatial frequency of -D^sd= ^ spatially Nyquist sampled 45 m telescope 

map. Therefore, the usage of the 45 m maps works to recover the flux densities missed by 
the mm-array. The basics of the data combination are described in Vogel et al. (1984), and 
our method here is similar to those in the previous works (e.g., Takakuwa et al. 2003). We 
took the following steps to make the combined maps. 
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1. For the BEARS H^^CO+ data, we used the spectra taken with the inner 9 beams cen- 
tered on the central beam on axis. The ahgnment of the 9 beams can be approximated 
as a regular grid with 41'/4 spacing for our RMS pointing accuracy of ^ 2". This makes 
the resultant map grid spacing of 20'.'7, i.e., full-beam sampling. Obtaining the data on 
a regular grid is required to use Fast Fourier Transformation (FFT). As for the N2H"'", 
we used the data on a regular grid of 17" spacing (full-beam sampling) with pointing 
accuracy of <^ 2". To achieve the excellent pointing accuracy, we used the spectra 
under wind speeds of ^ 2 m s~^. We corrected the scale of the H^^CO+ data with 
the BEARS in double sideband mode by observing the core center with the S80/S100 
receivers having single sideband filters. Since Tmb at the core center is 0.73±0.04 K 
by the S80/S100 receivers with 18'.'5 beam and 0.91±0.1 K by the BEARS with 2Qf!3 
beam, we scaled the intensities of the BEARS data by a factor of 0.80. 

2. To have a common frequency resolution between the OVRO and 45 m telescope data, 
the 45 m N2H+ spectra, which have a 37 kHz spectral resolution, are smoothed in the 
same frequency bins as those of the OVRO data with a 125 kHz resolution. For the 
H^^CO"*" hne, the 45 m and OVRO data have the same frequency bins. 

3. We made the 45 m maps by full-beam sampling on the regular grids with spacings of 
20'.'7 (= 41'/4/2) for R^^CO+ and 17'!0 for N2H+. We did not apply the task BPOSC 
in the NEWSTAR package, which corrects the misalignment between the 25 beams on 
the sky plane, to the H^3co+ data. Note that our OVRO field of view (FoV) with a 
single pointing does not cover the whole extent of the core. 

4. In converting the Tmb scale of the 45 m maps into the flux density per beam solid 
angle (F^,), we calculated conversion factors through a relation of Fj, = lu^heam 
with the Rayleigh- Jeans approximation. Here Ii, is the specific intensity and f^beam 
is the beam sohd angle. The conversion factors are given by ^ = || • ^^^^^^ — 



5. To deconvolve the 45 m maps with its beam, the beam pattern is approximated by 
an axisymmetric Gaussian profile with HPBW of 20'.'3 at 86.754 GHz for the BEARS 
H13C0+ data and 17"2 at 93.173 GHz for the S80/S100 N2H+ data. We adopt a mean 
HPBW of the inner 9 beams for BEARS (20'.'3 ± 0'/3) }i^^CO+ observations. Note that 
the S80/S100 receivers have nearly equal beam pattern as they detect two independent 
linearly polarized waves by sharing common dewar and waveguide. 




6. The 45 m maps are deconvolved by the Gaussian beams defined in (5). Since our 45 m 
maps satisfy the criteria discussed in Vogel et al. (1984), the deconvolution is done 
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by dividing the maps by the beam pattern. After several trials including the steps 
described below, the most plausible combined maps are obtained when the cut-off 
level of deconvolution is set to be 20% of the RMS noise level of the single-dish map 
so as not to miss high spatial frequency components having low amplitudes. 

7. The deconvolved 45 m maps are convolved with the primary beam pattern of the OVRO 
10.4 m antennas, assuming that all the six element telescopes have an axisymmetric 
Gaussian profile with HPBWs of 69'.'7 at 93.174 GHz and 75'/0 at 86.754 GHz. 

8. The single-dish maps are converted into visibihty data in the u-v domain by FFT. As 
demonstrated in Figure 16, the visibility amplitudes from the two instruments agree 
well with each other over the common u-v distances. We stress that we did not apply 
amphtude scale corrections to both the 45 m and OVRO data. 

9. The FFTcd single-dish visibilities are used to make model visibilities which have the 
same integration time as the OVRO data (205 sec). Since our 45m maps are sampled 
with the full-beam spacing, ^g^, we made model visibilities between zero and Df^^^ 
— l/(2^gj5) on the basis of the samphng theorem. We calculated -D^sd— 6.10 and 5.08 
kX for the N2H"'" and H^^CO"*" data, respectively, and conservatively adopted smaller 
cut-off values of 5.3 and 4.5 kX to fill the central u-v holes in the OVRO data (Table 
3). 

10. We combined the 45 m and OVRO visibihties in the u-v plane, assuming that the 45 m 
model visibilities have a Gaussian distribution within a radius of Df^^^ in the u-v plane. 
We found that synthesized beams of the combined maps are controlled by the visibility 
weighting function, regardless of the distribution of single-dish model visibilities. 

11. To construct an image from the combined visibilities, we used the task IMAGR in 
the AIPS package where one can optimize the synthesized beam size and sensitivity 
of the resultant image by changing the robustness power of the weighting function. 

Figure 17 represents how the parameter affects the synthesized beam pattern and the 
dirty map (map affected by the side lobe pattern of the synthesize beam). Considering 
the trade-off between the beam size and the smoothness of the map, we selected the 
robustness power of —1 and CLEANed the dirty maps. We stopped the CLEAN loop 
before encountering the first negative CLEAN component because we are unable to 
find emission- free regions to estimate an RMS noise level, i.e., a threshold flux density 
to stop the loop, due to the extent of the core than the OVRO FoV. 
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B. Calculation of Column Density 

To estimate the total column density of molecular gas (A^tot) through observations of 
the rotational transition from the upper J + 1 to lower J level, we use the previously derived 
equation (e.g., Scoville ct al. 1986) which gives A^tot through the optical depth (r^,) and 
excitation temperature (Tex) of the probe molecule as follows. 



3A; (Tex + f) e^Pt 



/ hB.J(.J+l) 



Ntot = ^^Y^ ■ ^^^—^ ^ / T,dv (Bl) 

SttV^ J+l i_exp(-J^)^ 

where k denotes the Boltzmann constant, h the Planck constant, fi the permanent dipole 
moment of the molecule, and B rotational constant, respectively. We should keep in mind 
that, to derive Eq.(Bl), one has approximated the rotational partition function Q{Tf,^) as, 

QiTe.)^Aln.+ ^]. (B2) 



hB \ 3k 



B.l. N2H+ 

To apply Eq.(Bl) to the N2H+ line data, we have to find Tex and r,^. In the following, we 
describe how they are obtained from the N2H+ spectra through the hyperfine structure (HPS) 
analysis. In local thermodynamic equilibrium (LTE), an observed brightness temperature 
(Tb) is expressed by the following radiative transfer equation. 



7b = / • [JM - ^.(7bg)]{l - e-"-^}, (B3) 

where Tb is replaced by the brightness temperatures averaged over a main beam (Tmb) 
for single-dish observations, and over a synthesized beam (Tsb) for interferometric obser- 
vations. Here / denotes a beam filling factor, and Ju{T) the equivalent Rayleigh- Jeans 
temperature given by JviT) = ^ |exp — l} . Prom Eq.(B3), we fit the function of 
y{v) = C [1 — exp {—t{v)}] with a constant C to each N2H+ hyperfine (HP) component. 
The optical depth is given by. 



r{v)^ Ttot Si exp 

i=l 



2a? 



(B4) 
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as a function of velocity. Here Ttot denotes the total optical depth, i.e., a sum of the 
optical depths at the centers of the HF components, Sj the normalized relative intensity of 
the ith HF component taken from Cazzoli et al. (1985), 6vi the velocity offset of the ith 
component, and af = ^^"'§7™'''' with the velocity width in FWHM of the ith component 
(AfFWHM.i)- We assumed that the line broadening due to the thermal and non-thermal gas 
motions equally work for all the HF components since they are emanated from the same 
molecular ion in the same region. We thus simplified Eq. (B4) by adopting ai ~ a and 
At»FWHM,i = Ai'FWHM- Therefore, the HFS analysis gives Ttot, the center velocity {Vsys + 5vi), 
the velocity width (Ai>fwhm), and the constant C. Since C = /{J,^(Tex) — Ji/(Tbg) }, one 
can calculate Tex by. 



k 



(B5) 



Here we assumed / = 1.0 because the region of our interest is larger than the beam 
sizes. We now have the estimates of Tex, Ttot, and Ai>int at each observing point (each pixel 
for the combined data). Since the integral in Eq.(Bl) becomes 




71" 

TtotAVint, (B6) 



we have the desired expression to calculate the column density of N2H+ (A^N2H+) as follows, 
AT.^H. = 3.30 X 10 ^^^zi;^ [y^^) i.q ^ 3-1 )• (^7) 



Here we used = 3.4 D (Havenith et al. 1990) and B = 46586.702 MHz (Caselh, Myers 
& Thaddeus 1995). 



B.2. H^^CO+ 

We estimated the column density of H^^CO^ (-^hi3co+) with Eq.(Bl) giving the mean 
Tex and r^^g over the full extent of the core (Table 5) where Tj^q is a mean optical depth at 
the line center frequency of z/q. Since the radiative transfer equation of Eq.(B3) is written 
as Tb = / ■ [J,^(Tex) — Ju{Thg)] ■ ''K the integral of Ti, can be expressed as, 

^^y)dv:^^^ /^b(^)^;^ forr.~0. (B8) 

^ ^ l-e---o J,(Tex) - J.(Tbg) ^ ^ 
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Recall that H^3C0+ (1-0) in GF9-2 is optically thin (§4.1.2) and we assumed / = 1 
in §B.l. Combining Eqs.(Bl) and (B8), we obtain an expression to estimate the H^^CO+ 
column density from single-dish data as, 



N _ooo ,,inll (^ex + 0-69) 1 ( jTray^dv \ 

Here we took // = 4.07 D (Haese & Woods 1979) and B = 43377.17 MHz. To have a 
complementary expression for interferometric data, we converted the flux density (Fi,) per 
synthesized beam into a brightness temperature per beam (Tsb) through Tgb = ■ F„ 
where is a beam solid angle given by (^maj x ^min)- Therefore, we have 



N -P77,,ini3 (^e. + 0-69) 1 (Jh^Y ( I F^dv \ 

iVH.3co.- 3.77X10 i_,-4.WT.-j.(Te.)-J.(Tb,)-ll^J ll.OJykms-j ^'"^ 

(BIO) 



B.3. CCS 

Column density analysis of the CCS 43 — 32 line is similar with that for the H^'^CO+ line 
in the sense that we have only the representative values of Tex and over the entire core 
(Table 5). The observed CCS brightness temperature are converted into the column density 
through 



^ocs = ^.l^^l^ r3^- J.(rex)-J.(r.j i l.0Kkms-- ^ 



(Bll) 



Here we used = 2.81 D (Murakami 1990), and B = 6477.75036 MHz (Yamamoto 
et al. 1990). Although the rotational energy [E-e) of the CCS molecule is given by Ej- = 
hBJ{J + 1) — hDJ'^{J + 1)^ considering the centrifugal stretching with D = 1727.96 kHz 
(Yamamoto et al. 1990), we ignored the D term to use the rotational partition function of 
Eq.(B2). 



-33- 



REFERENCES 

Aikawa, Y., Ohashi, N., Inutsuka, S., Herbst, E., & Takakuwa, S. 2001, ApJ, 552, 639 
Andre, P., Martin-Pintado, J., Despois, D., & Montmerle, T. 1990, A&A, 236, 180 
Andre, P., Ward-Thompson, D., & Barsony, M. 1993, ApJ, 406, 122 (AWB93) 
Andre, P., Ward-Thompson, D., & Motte, F. 1996, A&A, 314, 625 
Andre, P., Motte, F., & Bacmann, A. 1999, ApJ, 513, L57 

Andre, P., Ward-Thompson, D., & Barsony, M. 2000, Protostars and Planets IV, 59 
(AWBOO) 

Bachiller, R. 1996, ARA&A, 34, 111 

Beckwith, S. V. W., Henning, T., & Nakagawa, Y. 2000, Protostars and Planets IV, 533 

Belloche, A., & Andre, P. 2004, A&A, 419, L35 

Benson, P. J., CaseUi, P., & Myers, P. C. 1998, ApJ, 506, 743 

Bergin, E. A. & Langer, W. D. 1997, ApJ, 486, 316 

Butner, H. M., Lada, E. A., & Loren, R. B. 1995, ApJ, 448, 207 

Cabrit, S., & Bertout, C. 1992, A&A, 261, 274 

CaseUi, P., Myers, P. C, & Thaddeus, P. 1995, ApJ, 455, L77 

CaseUi, P., Benson, P. J., Myers, P. C, & Tafalla, M. 2002, ApJ, 572, 238 

Claussen, M. J., Marvel, K. B., Wootten, A., & Wilking, B. A. 1998, ApJ, 507, L79 

Chen, H., Myers, P. C, Ladd, E. F., & Wood, D. O. S. 1995, ApJ, 445, 377 

Chernin, L. M. 1995, ApJ, 440, L97 

Codella, C, Bachiller, R., & Reipurth, B. 1999, A&A, 343, 585 
Cox, P., Walmsley C. M., & Gusten, R. 1989, A&A, 209, 382 
Dickman, R. L. 1978, ApJS, 37, 407 
DoweU, C. D., et al. 2003, Proc. SPIE, 4855, 73 



Foster, P. N., & Chevalier, R. A. 1993, ApJ, 416, 303 

Furuya, R. S., Kitamura, Y., Saito, M., Kawabe, R., & Wootten, H. A. 1999, ApJ, 525, 821 

Furuya, R. S., Kitamura, Y., Wootten, H. A., Claussen, M. J., Saito, M., Marvel, K. B., & 
Kawabe, R. 2000, ApJ, 542, L135 

Furuya, R. S., Kitamura, Y., Wootten, H. A., Claussen, M. J., & Kawabe, R. 2001, ApJ, 
559, L143 

Furuya, R. S., Kitamura, Y., Wootten, A., Claussen, M. J., & Kawabe, R. 2003, ApJS, 144, 
71 

Goodman, A. A., Benson, R J., Fuller, G. A., & Myers, R C. 1993, ApJ, 406, 528 
Haese, N. N., & Woods, R. C. 1979, J. Chem. Phys., L61, 396 

Havenith, M., Zwart, E., Leo Meerts, W., & Ter Meulen, J. J. 1990, J. Chem. Phys., 93, 
8446 

Ho, P. T. P. & Townes, C. H. 1983, ARA&A, 21, 239 
Hunter, C. 1977, ApJ, 218, 834 

Jijina, J., Myers, P. C, & Adams, F. C. 1999, ApJS, 125, 161 

Kovacs, A. 2006, SHARC-2 350 //m Observations of Distant Submillimeter Galaxies and 
Techniques for the Optimal Analysis and Observing Weak Signals, Ph.D. Dissertation, 
Cahfornia Institute of Technology 

Lai, S.-R, & Crutcher, R. M. 2000, ApJS, 128, 271 

Larson, R. B. 1969, MNRAS, 145, 271 

Leong, M. M. 2005, URSI Conf. Sec. J3-10, 426 

Mizuno, A., Onishi, T., Hayashi, M., Ohashi, N., Sunada, K., Hasegawa, T., & Fukui, Y. 
1994, Nature, 368, 719 

Moriarty-Schieven, G. H., Wannier, P. G., Keene, J., & Tamura, M. 1994, ApJ, 436, 800 

Moscadelh, L., Testi, L., Furuya, R. S., Goddi, C, Claussen, M., Kitamura, Y., & Wootten, 
A. 2006, A&A, 446, 985 

Mundy, L. G., Looney, L. W., & Welch, W. J. 2000, Protostars and Planets IV, 355 



-35- 

Murakami, A. 1990, ApJ, 357, 288 

Myers, P. C, & Benson, P. J. 1983, ApJ, 266, 309 

Myers, P. C, & Ladd, E. F. 1993, ApJ, 413, L47 

Ogino, S., Tomisaka, K., & Nakamura, F. 1999, PASJ, 51, 637 

Ossenkopf, V., & Henning, T. 1994, A&A, 291, 943 

Penston, M. V. 1969, MNRAS, 144, 425 

Phillips, T. G., Wannier, P. G., Scoville, N. Z., & Huggins, P. J. 1979, ApJ, 231, 720 

Preibisch, T., Ossenkopf, V., Yorke, H. W., & Henning, T. 1993, A&A, 279, 577 

Richer, J. S., Padman, R., Ward- Thompson, D., Hills, R. E., & Harris, A. 1. 1993, MNRAS, 
262, 839 

Richer, J. S., Shepherd, D. S., Cabrit, S., Bachiller, R., & Churchwell, E. 2000, Protostars 
and Planets IV, 867 

Schilke, P., Walmsley, C. M., Pineau des Forets, G., & Flower, D. R. 1997, A&A, 321, 293 
Schneider, S., & Elmegreen, B. G. 1979, ApJS, 41, 87 

Scoville, N. Z., Sargent, A. I., Sanders, D. B., Claussen, M. J., Masson, C. R., Lo, K. Y., & 
Phillips, T. G. 1986, ApJ, 303, 416 

Shinnaga, H., Ohashi, N., Lee, S., & Moriarty-Schieven, G. H. 2004, ApJ, 601, 962 

Shu, F. H. 1977, ApJ, 214, 488 

Shu, F. H., Adams, F. C., & Lizano, S. 1987, ARA&A, 25, 23 
Stutzki, J., Ungerechts, H., & Winnewisser, G. 1982, A&A, 111, 201 

Suzuki, H., Yamamoto, S., Ohishi, M., Kaifu, N., Ishikawa, S., Hirahara, Y., & Takano, S. 
1992, ApJ, 392, 551 

Tafalla, M., Myers, P. C., Caselh, P., Walmsley, C. M., & Comito, C. 2002, ApJ, 569, 815 
Tafalla, M., Myers, P. C., Caselh, P., & Walmsley, C. M. 2004, A&A, 416, 191 
Terebey, S., Vogel, S. N., & Myers, P. C. 1992, ApJ, 390, 181 



-36- 

Takakuwa, S., Kamazaki, T., Saito, M., & Hirano, N. 2003, ApJ, 584, 818 
Ungerechts, H., Walmsley, C. M., & Winnewisser, G. 1980, A&A, 88, 259 
Vogel, S. N., Wright, M. C. H., Plambeck, R. L., & Welch, W. J. 1984, ApJ, 283, 655 
Yamamoto, S., ct al., 1990, ApJ, 361, 318 

Walker, C. K., Lada, C. J., Young, E. T., Maloney, P. R., & Wilking, B. A. 1986, ApJ, 309, 
L47 

Ward-Thompson, D., Scott, R F., Hills, R. E., & Andre, P. 1994, MNRAS, 268, 276 
Whitworth, A., & Summers, D. 1985, MNRAS, 214, 1 

Wiesemeyer, H. 1997, The spectral signature of accretion in low-mass protostars, Ph.D. 
Dissertation, University of Bonn 

Wiesemeyer, H., Cox, P., Giisten, R., & Zylka, R. 1999, ESA SP-427: The Universe as Seen 
by ISO, 533 

Winnewisser, G., Churchwell, E., & Walmsley, C. M. 1979, A&A, 72, 215 
Wolkovitch, D., Langer, W. D., Goldsmith, P. P., & Heyer, M. 1997, ApJ, 477, 241 
Womack, M., Ziurys, L. M., & Wyckoff, S. 1992, ApJ, 387, 417 
Zhou, S. 1992, ApJ, 394, 204 



This preprint was prepared with the A AS IM^jX macros v5.2. 



-37- 



Table 1. Summary of Single- Dish Telescope Observations 



Emission 


Rest Frequency 
(MHz) 


Xelescope 


Receiver 


^HPBW* 

(circsec) 


'?mb 


(km s-l) 


(mK) 


Type 


Area^ 
(arcmin) 


NH3 (1,1) 


23694.495 


NRO 


H22 


78 


0.82 


0.477 


46 


M/N 


5.9 X 5.9 


NH3 (2,2) 


23722.633 


NRO 


H22 


78 


0.82 


0.476 


56 


M/N 


5.9 X 5.9 


NH3 (3,3) 


23870.130 


NRO 


H22 


78 


0.82 


0.473 


58 


M/N 


5.9 X 5.9 


NH3 (4,4) 


24139.417 


NRO 


H22 


78 


0.82 


0.469 


62 


M/N 


5.9 X 5.9 


CC^'^S 43 - 39 


44497.559 


NRO 


840 


40 


0.76 


0.253 


25 


D 




CCS 43 - 32 


45379.033 


NRO 


840 


40 


0.76 


0.245 


50 


M/F 


5.3 X 5.3 


HCi*0+ (1-0) 


85162.157 


NRO 


880 


18.7 


0.50 


0.130 


26 


D 




C3H2 2i2 — l()l 


85338.906 


NRO 


8100 


18.7 


0.50 


0.130 


30 


D 




H13CO+ (1-0) 


86754.330 


NRO 


S80/S100 


18.5 


0.50 


0.128 


37 


C/F 


3.7 X 2.4 




86754.330 


NRO 


BEARS 


20.3 


0.49 


0.108 


48 


M/F 


8.6 X 8.6 


SiO (1-0) v = 


86846.998 


NRO 


880 


18.5 


0.50 


0.128 


42 


M/F 


3.7 X 2.4 




86846.998 


NRO 


BEARS 


20.5 


0.49 


0.108 


185 


M/F 


8.6 X 8.6 


N2H+(l-0) 


93173.900 


NRO 


S80/S100 


17.2 


0.51 


0.119 


31 


M/F 


2.0 X 2.0 


12 CO (1-0) 


115271.204 


NRO 


8100 


14.8 


0.49 


0.0976 


94 


C/F 


2.8 X 1.8 




115271.204 


NRO 


BEARS 


16.5 


0.51 


0.0983 


139 


M/F 


7.1 X 7.1 


12 CO (3-2) 


345795.989 


CSO 




22.0 


0.75 


0.416 


290 


M/F 


3.7 X 3.0 


350 /nm continuum 


850000 


CSO 


SHARCII 


8.9x8.0 










4.0 X 3.8 



^Half-power beam width for a circular Gaussian beam 
''Main beam efficiency 
'^Effective velocity resolution 

"^Typical RMS noise level of the spectrum. For detection observations, an RMS noise level of the spectrum is given. For mapping 

observations, a mean RMS noise level for all the spectra is given. 

■^Observing modes for the molecular line emission; D denotes "detection" observations toward the center position of R.A.(J2000)= 
20''51'"29.5'', Decl.(J2000)= 60'*18'38.0", M/F and M/N denote, respectively, "mapping" observations with full-beam and Nyquist 
sampled grid spacings centered on the above position, and C/F denotes "cross scan" observations with a full-beam spacing along the 
NE-SW and NW-SE directions centered on the position. 

^Size of the region for the "mapping" observations, and the major and minor axis lengths for "cross scan" ones. 



Table 2. Summary of Interferometric Observations 







Center 


Synthesized Beam 


Bandwidth 


Velocity 


Image 


Array 


Emission 


Frequency 


^maj ^ ^min 


RA. 


Resolution 


Sensitivity 






(MHz) 


(arcsecxarcsec) 


(deg) 


(MHz) 


(km s-i) 


(mjy beam"') 


VLA 


continuum 


8435 


3.29x2.93 


29.7 


172* 




0.038 


OVRO 


continuum 


89964 


5.38x4.80 


41.3 


4096'' 




0.15 




N2H+ (1 0) 


93173.78 


4.91x4. .36 


-62.4 


7.75 


0.402 


17,6^ 




Il^-'CO 1 (1 0) 


,S(i751. :-!:-! 


().:-!! X 1.81 


-21. V 


l,!):-!8 


0,108 


:i--!,0' 



''Concatenated the two polarization data. 
''Concatenated the lower and upper sideband data. 
"Per velocity channel. 
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Table 3. Summary of the Combined Data from the OVRO MM- Array and the NRO 45 M 

Telescope^ 



Line 



Spatial Frequency Range 



45 m 
(kA) 



OVRO 
(kA) 



Field of 

View 
(circsec) 



Synthesized Beam 

maj ^ ^min P.A. 

(arcsec) (deg) 



Velocity 
Resolution 
(km s-i) 



N2H+(l-0) 0.0-5.29 3.30-71.8 

H"CO+ (1-0) 0.0-4.48 3.07-58.0 



69.7 
75.0 



4.49x3.86 
6.73x5.14 



-70.7 
-27.2 



0.402 
0.108 



"For details, see Appendix A 



Tabic 4. Results of Contimmm Emission Fhix Measurements 







Size 




Peak 


Total Flux 


A 


Instruments 




P.A. 


Intensity 


Density 


(/im) 




(arcsec x arcsec) 


(deg) 


(mjy beam~^) 


(mJy) 


35540 


VLA-D 






< 0.055 




3330 


OVRO 


(6.1 ±0.9) X (4.3 ±0.6) 


129±17 


1.16±0.17 


0.75±0.07 


350 


CSO-SHARCII 


(32.9 ±4.3) X (22.5 ±2.9) 


140±13 


300±30 


1830±330 
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Table 5. Results from the 45 m Telescope Spectral Line Observations: Dense Gas Tracers 



Core Center* Inside the Half-Maximum Contour'' 



Line 
















(K) 


(km s^^) 




(km s-l) 




(K) 


N2H+(l-0) 


1.96±0.05 


0.35±0.01 


3.4±0.08 


0.32±0.09 


2.7±0.4 


9.5±1.9 


H13CO+ (1-0) 


0.88±0.09 


0.60±0.03 


< 0.1 


0.60± 0.17 






C3H2 2l2 — loi 


1.12±0.06 


0.59±0.1 


0.34 








CCS 43 - 32 


0.74±0.07 


0.41±0.04 


0.94 


0.41±0.19 






NHa (1,1) 


1.20±0.06 


0.61±0.10 


1.6±0.1 


0.62±0.17 


0.45±0.18 


< 9.7 



''Single pointing toward the 3 mm continuum source. See the spectra in Figure 2, and beam sizes 
in Table 1. 

''For N2H+ and NH3 lines which allow to derive optical depth (r) and excitation temperature 
(Tex) through hyperfine structure (HFS) analysis, we present mean and standard deviation for the 
spectra existing inside the half-maximum, i.e., 50% level contour of total integrated intensity map 
(Figure 4). 

'^Peak main-beam brightness temperature. For N2H+ and NH3 lines, those for the brightest HF 
components are given. 

''in calculating intrinsic velocity width (Avi^t), we employed HFS analysis for the N2H+ and 
NH3 lines, spectral moment analysis for the H'^CO+ and C3H2 lines, and single Gaussian fitting 
for the CCS line (§4.5.2). For the C3H2 and CCS lines, we corrected for line broadening due to the 
optical depth (Phillips et al. 1969). 

"Optical depth. We present t of the brightest hyperfine transition Fi,F = 2,3—1,2 for N2H+ 
(1-0) (§4.1.1) and F,Fi= f ,2- f,2 for NH3 (1,1) (§4.1.5). 
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Table 6. Column Densities of N2H+, Hi3co+, C3H2, CCS, and NH3 Molecules^ 





Combined 




Nobeyama 45 m Telescope 




Line 


Inside 350 /xm 




Core Center"^ 




Pull Extent of 




3a Contour'' 








the Core'' 






r < 10" 


r < 20" 


r < 40" 






(cm- 2) 


(cm- 2) 


(cm-2) 


(cm-2) 


(cm-2) 


N2H+ 


(5.2±2.5)E+13 


(5.6±0.2)E+13 


(4.4±0.3)E+13 


(3.5±0.2)E+13 


(3.9±1.0)E+13 


Hi-^CO+ 


(4.7±0.8)E+11 


(6.1±0.3)E+11 






(4.3±0.2)E+11 


C3H2 




(2.4±0.3)E+13 








CCS 






(1.5±0.2)E+13 




(1.5±0.3)E+13 


NH3 








(8.2±0.2)E+14 


(7.7±2.8)E+14 



''Calculated with the assumption of Tex = 9.5 K (§4.1.1) 

''Column densities measured inside an area enclosed by the 3<T level contour of the 350 fini continuum emission. 
For this purpose, we made the N2H+ and H'^CO"'" maps by giving the SHARC II beam size as a synthesized 
beam size. 

"^Toward the 3mm continuum source measured witli the beam size in Table 1. e.g., 9hpbw = 20" beam for 
the region of r < 10". Regarding the N2H+ column densities for the regions of r < 20" and 40", we calculate 
column densities using these spectra smoothed with ^hpbw = 40" and 80" beams, respectively (§4.1.1). 

''Measured over the full extent of the core with beam sizes of 20" for N2H+, H'3C0+, and C3H2, 40" for 
CCS, and 80" for NH3. Here the full extent of the core is defined by a region enclosed by the 50% level contour 
of each total intensity map. 



Table 7. Fractional Abundances of Dense Gas Tracer Molecules 



Beam Size 

Molecule Xm* (arcsec) Comments 



N2H+ (9.5 ± 5.0)E-10 9x8 Used {N(H2))'° in §3.1.2 

H13CO+ (8.5 ± 4.7)E-12 9x8 Used (Ar(H2))'' in §3.1.2 

C3H2 (2±1)E-10 20 X(N2H+)^„ X fJ^^(20")'^ 

CCS (5±3)E-10 40 X(N2H+)^„ X ^1^1^(40")- 

NH3 (2±l)E-8 80 x(N2H+)_„ X (80")= 



''■Fractional abundance of the molecule. 

''We adopted mean H2 column density, (iV(H2)), from cases I and II in §3.1.2 
and considered the range of H2 column densities into the errors in Xm. 

"^Assumed that N2H+ does not have radial abundance variation. Column 
column density ratios are calculated from Table 6. The notation of X(N2H+)g// 
indicates N2H+ abundance measured with the SHARC II beam. 
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Table 8. LTE and Virial Masses of GF 9-2 Core 















Line 






Res" 


Case l" 


Case IP 




(pc^) 




(pc) 






N2H+(l-0) 


2.632E-3 


1.2±0.7 


0.0289 


0.39±0.3 


O.lOiO.08 


H13CO+ (1-0) 


3.963E-3 


1.3±0.3 


0.0354 


2.9±1.7 


0.12±0.07 


CCS 43 - 32 


5.897E-3 


3.1±1.9 


0.0432 


1.7±1.4 


0.078±0.06 


NHs (1,1) 


5.544E-3 


4.5±2.4 


0.0419 


2.3±1.6 


0.25±0.15 



''Area enclosed by the 50% level contour of each total integrated intensity map 
(Figure 4). 

''Errors in LTE-mass include those in the column density, area size, and the 
abundance. 

'^Effective radius defined by i?eff = \/ A/tt. 

''Errors in virial mass are calculated from (<5Afvir)^ = ^ d(Av"t} ) (^("^^int))^ + 

( a Jt^f f ) ("^^efl)^ where 6i?gff and (5(At!int) are those in R^g and (Atiint), respec- 
tively. The intrinsic velocity width (Avint) and its error are taken from those for 

the entire core in Tabic 5. 

"Cases I and II are, respectively, when (Aiiint)^ = (Atithm)^ + ('^^turb)'^ ^^'^ 
(Ai^int)^ = (A^Jthm)^- Here (Ai^int), (At)thm>, and (Avturb) are means of intrinsic, 
thermal, and turbulent velocity widths, respectively (see §5.2.1). 
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Fig. 1. — Overlay of the 3 mm continuum emission map observed with the OVRO mm- 
array (grey scale plus contours) on the 350 //m map with the CSO SHARC II camera (thick 
contours). The contours start from the 3a levels with a 3a interval where cr = 0.15 and 29.7 
mJy beam~^ for the 3 mm and 350 //m images, respectively. The dashed contours show the 
—3a level for the 3 mm image. The 350 /xm image docs not have negative contours below the 
—2.2a level. The ellipses at the bottom corners indicate the beam sizes. 



Fig. 2. — Spectra of the molecular line emission shown in the main-beam brightness temper- 
ature (Tmb) scale toward the 3 mm source. The vertical dashed line indicates the systemic 
velocity (Kys) of —2.5 km s~^ (§3.2.1). The horizontal blue and red bars under the ^^CO 
lines indicate the LSR-velocity ranges for the integrated intensity maps shown in Figure 8, 
and the green bars under the N2H+ and H^^CO"*" lines are for the maps shown in Figure 6. 
The N2H+ emission seen in 4.5 ^ li,sR/km s~^ ^ 7.2 is the main hyperfine components (see 
Figure 3). 



Fig. 3. — N2H"'" (1-0) spectrum (black histogram) taken with the 45 m telescope toward the 
3 mm source. The green histogram represents the best-fit spectrum from the HF analysis 
(§4.1.1). We assume that the LSR velocity of the isolated Fi,F=0,l-l,2 transition (z^rcst = 
93176.310 MHz) is equal to the Vgys of —2.5 km s^^. The lengths of vertical lines above the 
spectrum indicate the expected relative intensities at r ~ 0, plotted at the rest frequencies 
(CazoUi et al. 1985). The horizontal bars below the spectrum show the velocity ranges 
integrated to obtain maps in Figure 4c and those in the first row of Figure 6. 
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Fig. 4.— Total integrated intensity maps of the (a) Hi3C0+ (1-0), (b) CCS 43-82, (c) N2H+ 
(1-0), and (d) NH3 (1,1) lines observed with the 45m telescope. The velocity ranges for the 
integrations are between Vlsr = —2.95 and —2.1 km s~^ for H^^CO+, and —3.0 and —2.1 
km for CCS. For the N2H+ and NH3 emission, all the hyperfine emission are integrated 
over the velocity ranges shown by the horizontal bars in Figures 3 and 10, respectively. 
Contour intervals are 3(7, starting from the 3a levels. The la noise levels are 33.4, 42.5, 27.0, 
and 24.8 mK-km s"^ for the H^^CO+, CCS, N2H+, and NH3 maps, respectively. The central 
star marks the position of the 3 mm source. The hatched circles at the bottom- left corner 
of each panel show beam size (see Table 1), and the small dots represent the observed grid 
points. 



Fig. 5.— Velocity channel maps of the N2H+ (1-0) Fi,F=0,l-l,2 and Hi3C0+ (1-0) emission 
taken with the 45 m telescope (see Table 1); each panel has a size of 225" x 225". The central 
LSR-velocity for each channel is shown at the upper right corner. Contour intervals are 2a, 
starting from the 3a levels of the images where la noise level is 78 mK in for the N2H+ 
and 96 mK for H^^CO+. The star in each panel is the same as in Figure 4. The systemic 
velocity is Vj^sr — —2.5 km s~^. 

Fig. 6. — Comparison of the total integrated intensity maps of all the seven hyperfine (HF) emission 
of N2H+ (1-0) (first row), N2H+ (1-0) Fi,F=0,l-l,2 (second row) emission, and Hi3C0+ (10) 
(third row) line taken with the Nobeyama 45m telescope (first column), the OVRO mm-array 
(second column), and made from the combined 45 m plus OVRO visibilities (third column). The 
central star and the ellipse at the bottom left corner in each panel are the same as in Figure 4. The 
large circles in the second and third columns indicate the field of view. The N2H+ total emission 
maps are integrated over the 3 velocity ranges shown by the horizontal bars in Figure 3, and those 
for N2H"'' Fi,F = 0,1-1,2 and H^'^CG"*" maps by the green ones in Figure 2. The contour intervals 
are the same as in Figure 1. The la rms noise levels are 38.5 mK in T^b for the 45 m N2H"'" total 
emission map, 42.5 mK for the 45 m N2H+ Fi,F = 0,1-1,2 map, 33.5 mK for the 45 m 11^^00+ 
map, 17.4 mJy beam-^ for the OVRO N2H+ total emission image, 22.4 mJy beam ^ for the OVRO 
N2H+ Fi,F = 0,1-1,2 image, 17.2 mJy beam-^ for the OVRO Hi3C0+ image, 62.4 mJy beam-^ 
for the combined N2H"'" total emission image, 95.5 mJy beam^^ for the combined N2H"'" Fi,F = 
0,1-1,2 image, and 35.3 mJy beam^^ for the combined H^'^CO^ image. The noise levels for the 
combined data correspond to 170 mK in for the H^^CO"*", 510 mK for the N2H"'" total, and 780 
mK for the N2H''" Fi,F = 0,1-1,2 emission maps. 
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Fig. 7. — Velocity channel maps of the combined OVRO plus 45 m telescope images of the 
N2H+ (1-0) Fi,F=0,l-L2 and Hi-^CO+ (1-0) hnes (see Table 3). Contour intervals are the 
same as in Figure 1, and the RMS noise levels are 32.1 and 52.1 mJy beam^^ for the N2H^ 
and H^^CO+ lines, respectively. The size for each channel map is 76'/8x76'/8. The systemic 
velocity is Vlsr = —2.5 km s""^. 



Fig. 8. — Maps of the blue- (left panels) and redshifted (right panels) -"^^CO (1-0) (upper panels) 
and (3-2) (lower panels) emission. The (1-0) emission is integrated over the velocity ranges of 
-4.2 < T^LSR/km < -3.4 for the blue and -1.8 < lisR/km s"^ < -1.0 for the red. The 
(3-2) emission is integrated in the ranges of —5.8 < ^LSR/km s^^ < —3.8 for the blue and —1.4 < 
^LSR/km s~^ < +0.6 for the red. These velocity ranges are shown in Figure 2. Contours with 
color scale are drawn by the 10% step up to the 90% level of the peak value. The lowest contour 
for the (1-0) maps is the 10% level, corresponding to the 21a and 8.4cr levels for the blue and red, 
respectively. In the (3-2) maps, the lowest contour is the 50% level, corresponding to the 9. Oct and 
8.80" levels for the blue and red, respectively. The white dashed contours shown in all the panels 
indicate the 5cr level contour of the H^'^CO^ total intensity map (Figure 4a), and the inner ones 
in (c) and (d) indicate the 50% (corresponding to 16.40") level contour. The central circles with a 
radius of 40" in (a) and (b) show the error circle of the H2O maser position (Puruya et al. 2003). 
The central yellow star indicates the position of the 3 mm source. The other symbols are the same 
as in Figure 4. 



Fig. 9. — Column density maps of the N2H"'" line obtained from (a) the 45 m telescope data 
and (b) the combined OVRO plus 45 m data (§4.1) in color scale. The color scale bars are 
the same for both the 45 m and combined data. The black contours in (a) and white ones 
in (b) indicate the total integrated intensity maps of the N2H+ emission (Figure 4c) and the 
350 /um continuum emission (Figure 1), respectively. The white circle in (a) and black one 
in (b) indicate the OVRO filed of view. 
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Fig. 10. — NH3 (1,1) and (2,2) spectra observed with the 45 m telescope toward the 3 mm 
continuum emission. In the lower panel, the green histogram represents the best-fit NH3 (1,1) 
spectrum from the HFS analysis (§4.1.5). It is assumed that the (1,1) and (2,2) spectra have 
their brightest HP components at Flsr — —2.5 km s""^. The vertical lines above the spectra 
indicate the expected relative intensities at r ~ and LSR-velocities of the HF components. 
The horizontal bars below the (1,1) spectrum show the velocity ranges to obtain the total 
intensity map shown in Figure 4d. 



Fig. 11. — Radial column density profiles of the N2II+ in GF9-2 core obtained from the 
Nobeyama 45 m telescope (blue) and the combined OVRO plus 45 m data (black) (see text). 
The upper horizontal axis shows the linear scale. The right hand vertical axis represents 
H2 column density converted with X(N2lI+) = 9.5E-10 (Table 7). The solid curves indicate 
the best-fit power-law profiles (§4.4): N^,n+{r) = (7.3 ± 1.1) x lO^^ (r/l'.'0)"^^-°^°-^^ cm-^ for 
r < 50" for the combined and (8.5 ±2.5) x 10^^ (r/l'.'0)~^^-^^°-^^ cm^^ for r < 80" for the 45 m 
data. The broken curves indicate the beam patterns: 17"2 in HPBW for the 45 m telescope, 
and 4"1 for the combined data. The vertical dotted lines at r = 29'.'9 (corresponding to 0.029 
pc) indicates the Res for the N2H+ core (Table 8). 
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Fig. 12.— Position- Velocity (PV) diagrams of the N2H+ (1-0) Fi,F=0,l-l,2 (upper panels) 
and H^^CO+ (1-0) (lower panels) emission for GF9-2 core made from the 45 m telescope 
(green contour) and the combined data (greyscale with black contour). The left and right 
panels, respectively, show the diagrams along the major (P.A.= 45°) and minor (P.A.= 135°) 
axes, passing through the continuum source. The right hand vertical axis in each panel shows 
distance from the 3 mm source in pc. The contour intervals are the same as in Figure 4. The 
intensity is presented in T^h scale for the 45 m and in for the combined data. See the 
caption of Figure 6 for the RMS noise levels. The vertical and horizontal dotted lines indicate 
V^ys = —2.5 km s~^ and the position of the 3 mm source, respectively. The horizontal dashed 
lines in (c) and (d) running at ±37'/5 indicate the field of view of the combined data. The 
horizontal and vertical bars labeled with Vres and ^beam indicate the velocity resolution and 
beam size, respectively. 



Fig. 13. — Intrinsic velocity width (Avi^t) maps (color scale) of (a) N2H"'" (1-0) Fi,F=0,l-l,2 
and (b) H^^CO"'" (1-0) emission for the 45m telescope data (§4.5.2). The presented Avi^^t 
ranges are 0.28 < Afjnt/km s^^ < 0.61 for the N2H"'" and 0.44 < Atiint/km s~^ < 0.63 for the 
H^^CO^. Contours in (a) and (b) indicate the N2H"'" and H^^CO^ total integrated intensity 
maps shown in Figure 6. The thick white contour in each panel indicates the 50% level 
contour of the total maps. The star and hatched eUipse are the same as in Figure 4. 



Fig. 14. — Plot of normalized CO outflow momentum flux (Fco • c/Lboi) vs. bolometric 
temperature (Tboi) for nearby {d < 350 pc) low- mass YSOs less luminous than 5OL0 (§5.1), 
including S106FIR (Lboi ^ IO8OL0; Richer et al. 1993; Furuya et al. 2000). The fiUed- 
and open circles represent the presence and the absence of H2O masers in the period of the 
Nobeyama survey (Furuya et al. 2001; 2003), respectively. The locus of GF9-2 (Lboi — 
0.3 Lq) is indicated by large filled circle. The upper limit in Fqq " c/Lboi for S106FIR is 
recalculated from Furuya et al. (2000) with the 3a upper limit and the outflow velocity 
of 50 km s~^. The Tboi is taken from Chen et al. (1995) and AWBOO, and Fco ■ c/Lboi is 
calculated from Table 4 in Furuya et al. (2003). The Tboi ranges for the class 0, 1, and II 
sources (Chen et al. 1995) are shown by the bottom bars. 
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Fig. 15. — Map of the intrinsic velocity width of the H^^CO^ (1^0) emission (greyscale 
and contour), obtained from the combined data. Contour intervals are the 10% step of the 
difference between the maximum (0.71 km s~^) and minimum (0.36 km s~^) velocity widths. 
The thick white and black contours show the ISa level of the total integrated intensity map 
of the combined N2H+ data (Figure 6) and the 3a level of the 3 mm continuum emission 
(Figure 1), respectively. The green contours show the highest 3 level contours (18, 21 and 
24(7 levels) of the total integrated intensity map of the combined H^^CO"'" data (Figure 6). 
The ellipses are the same as in Figure 4. 



Fig. 16. — Upper and lower panels: plots of flux density in Jy vs. projected baseline length 
in kX for the N2H"'' (1-0) line data with the 45 m telescope (filled circles) and the OVRO 
mm-array (open triangles). The data are from Vlsr = —2.6 km s~^ channel of the strongest 
hyperfine component of I23 — O12 transition. The lower panel is a magnified plot of the upper 
one. The visibilities in the plots are resampled with 0.2 and 0.4 A;A bins for the 45 m and 
OVRO data, respectively. 



Fig. 17. — Beam patterns (upper panels) and "dirty maps" (lower panels) of the combined 
N2H+ (1-0) data for the robustness powers of +1, 0, —1, and —2. The data are from the —2.6 
km s~^ channel map of the strongest hyperfine component of I23 — O12 transition. Contour 
intervals for the beams are 10% of the peak values, starting from the ±10% levels. Dashed 
contours are complementary negative contours. Contours in the dirty maps have 3a steps, 
starting from the 3(7 levels; the la noise levels are 127.5, 82.4, 28.5, and 24.7 mJy beam~^ 
for the robustness powers of +1, 0, —1, and —2, respectively. The synthesized beam sizes 
of ^i„aj y<Omin are 10'.'33 x 9'.'18, 5'.'81 x 5'.'09, 4'.'16 x 3'.'57, and 3'.'27 x 2'.'66 for the robustness 
powers of +1, 0, —1, and —2, respectively. The central stars in the lower panels represent 
the position of the 3 mm source. 



This figure "fl.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



This figure "f2.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



This figure "f3.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



This figure "f4.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



This figure "f5a.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



This figure "f5b.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



This figure "f6.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



This figure "f7a.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



This figure "f7b.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



This figure "f8a.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



This figure "f8b.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



This figure "f8c.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



This figure "f8d.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



This figure "f9.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



This figure "flO.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



This figure "fll.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



This figure "fl2a.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



This figure "fl2b.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



This figure "fl2c.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



This figure "fl2d.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



This figure "fl3.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



This figure "fl4.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



This figure "fl5.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



This figure "fl6.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



This figure "fl7.jpg" is available in "jpg" format from: 



http ://arXiv. org/ps/astro-ph/0608 357 V 1 



